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II. On the Potential of Ellipsoidal Bodies, and the Figures of Equilibrium of 

Rotating Liquid Masses. 
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Eeceived May 29, — Eead June 25, 1914. 

By an ellipsoidal body is meant, in the present paper, any homogeneous body which 
can be arrived at by continuous distortion of an ellipsoid. If f = is the equation 
of the ellipsoid from which we start, and e is a parameter, the distortion of the 
ellipsoid may be supposed to proceed by e increasing from the value e = upwards, 
and the final figure may be taken to be 

fo+ef 1 + e% + e%+ ... - 0. 

For very small distortions the first two terms will adequately represent the 
distorted figure, and as we pass to higher orders the remaining terms will enter 
successively. 

The potential problem, to some extent interesting in itself, derives its chief 
importance from its application to the determination of the possible figures of 
equilibrium of a rotating mass of liquid. Poincare, # using his ingenious method 
of double layers, has shown how the potential of an ellipsoidal body can be carried 
as far as the second-order terms when the distortion is small, but gives no indication 
of how it is possible to carry it further, and indeed his method is one which hardly 
seems susceptible of being developed further than he himself has developed it. 
It is clear, however, that progress with the problem of rotating liquids can only 
be made when a method is available for writing down the potential of an ellipsoidal 
body distorted as far as we please. I believe the method explained in the present 
paper will be found capable of giving the potential of a body distorted to any 
extent, although (for reasons which will be explained later) I have not in the 
present paper carried the calculations further than second-order terms. 

The theory of figures of equilibrium of rotating masses of liquid is at present 
in an unsatisfactory state. It has been shown by Lord Kelvin that the Jacobian 
ellipsoid is stable at the point at which it coalesces with the Maclaurin spheroid, 
and it has been .shown by Poinoare to remain stable up to the point at which 

* " Sur la Stabilite de l'Equilibre des Figures Pyriformes affectees par une Masse Fluide en Kotation," 
< Phil Trans.,' A, vol. 198, p. 333. 

(524.) E 2 [Published February 2, 1915. 
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the series of Jacobian ellipsoids coalesces with the Poincare series of pear-shaped 
figures. After this point the series of Jacobian ellipsoids must, in accordance with 
Poinc are's doctrine of exchange of stabilities at a point of bifurcation, lose its 
stability, but the question of how it loses its stability is in a state of doubt. 
Darwin believed he had proved the Poincare series to be initially stable,^ whereas 
Liapounoff t has maintained that this series is initially unstable. The importance 
of this question to theories of cosmogony is, of course, great, although perhaps 
liable to be overrated. A caution of Poincare's J may be borne in mind : 
"Quelle que soit 1'hypothese [stability or instability] que doive triompher un jour, 
je tiens a mettre toute de suite en garde contre les consequences cosmogoniques 
qu'on pourrait en tirer. Les masses de la nature ne sont pas homogenes, et si on 
reconnaissait que les figures pyriformes sont instables, il pourrait neanmoins arriver 
qu une masse heterogene fut susceptible de prendre une forme d'equilibre analogue 
aux figures pyriformes, et qui serait stable. Le contraire pourrait d'ailleurs arriver 
egalement." 

The present investigation was started primarily in the hope of setting this 
question of stability at rest. I realised that to make a new series of computations 
on the subject could be of little value, for whatever the result, there would have 
been two opinions on the one side to one on the other. Moreover, Darwin has 
stated clearly that he does not think the divergence of opinion between 
M. Liapounoff and himself is one to be settled by new computations § : "I feel 
a conviction that the source of our disagreement will be found in some matter of 
principle/' I had hoped that it might be found possible to discuss the problem 
by a purely algebraical method, involving neither laborious computations nor 
intricate physical arguments, and that if such a discussion did not give a con- 
vincing and satisfying answer to the question in hand, at least it might reveal the 
source of disagreement between Darwin and Liapounoff. The result arrived at 
is one which, as will readily be understood from its nature, is only put forward 
with the utmost diffidence, but it is one from which I can find absolutely no escape. 
It is that underlying the whole question there is a complication, unsuspected equally 
by Poincare, Darwin, and (in so far as I can read his writings) Liapounoff, 
which renders nugatory the work of all these investigators on the stability of the 
pear-shaped figure. If my method is sound, it appears, as will be explained later, 
that it is impossible to draw any inference as to the stability of the pear from 
computations carried only as far as the second order of small quantities. The 

* " The Stability of the Pear-shaped Figure of Equilibrium of a Rotating Mass of Liquid," ' Phil. 
Trans./ 200 A (1902), p. 251; also papers in < Phil. Trans./ 208 A (1908), p. 1, and < Proc, Koy. Soc, 
82 A (1909), p. 188, all combined in one paper in 'ColL Scientific Papers/ vol. 3, p. 317. 

t "Sur un Probleme de Tchebychef," ' Memoires de TAcademie de St. Petersbourg/ xvii., 3 (1905). 

I Loc. cit. t p. 335. 

§ ' Coll. Scientific Papers/ 3, p. 392. 
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materials for an answer to the question are to be found only through the third- 
order terms. # Fortunately the method of the present paper admits of extension 
to the computation of third-order terms, and so it ought to, and I hope will be quite 
feasible to decide as to the stability or instability of the pear, a question reserved 
for a subsequent paper. 

The reader who is interested in the main conclusions of the paper rather than 
in details of theory, method, or calculations, may care to pass directly to § 35. 



General Theory of Potential op Ellipsoidal Bodies. 

2. We proceed to develop a method for writing down the potentials of certain 
homogeneous solids ; in particular of ellipsoids and distorted ellipsoids. We are for 
the present concerned solely with potential-theory — the discussion of rotating liquids 
does not enter before § 19. 

As will soon be evident, the problem in potential theory amounts to the following : 
to write the equation of the boundary of a homogeneous solid in such a form 
F (as, y, z) = 0, that the potential at the boundary is of the form F' (x, y, z) = 0, where 
F' (as, y, z) is a function containing the same algebraic terms as F (as, y, z), but having 
in general different coefficients. If this can be done, it only remains to equate 
F ; (as, y, z) + ^oo 2 (x 2 + y 2 ) to F (as, y, z\ and we have at once, on equating coefficients, 
a series of equations which will determine the possible figures of equilibrium for a liquid 
mass rotating with angular velocity «. 

3. Let F(as, y, z) = be the boundary of any homogeneous solid of density p. 
Assuming it to be possible, t let V t - be a function of position satisfying V 2 V,- = — iirp 
at all points of space and coinciding with the potential of the solid at all points inside 
the solid, and let V similarly be a function of position satisfying V 2 V = at all points 
of space, except possibly the origin or other infinitesimal region inside the solid, and 
coinciding with the potential of the solid at all points outside the solid. 

Then V t - must be equal to V at the boundary of the solid, and we must also have 

nr" " = T~° a * the boundary, where — denotes differentiation along the normal to the 

an an an 

surface. 

* Since writing this paper, I have been surprised to find that this conclusion is quite clearly implied 
in a paper which I published in 1902, "On the Equilibrium of Kotating Liquid Cylinders," 'Phil. 
Trans.,' A, 200, p. 67. See below, § 36. 

t I have not examined in any detail the conditions that this maybe possible, because the result of 
the paper proves that it is possible in the cases which are of importance. Similarly I ha^e not examined 
in detail the difficulties which might arise at the origin or at infinity, because in the final result they 
do not arise. We are searching for, and ultimately find, a certain solution of the potential equations, 
and after the solution has been obtained it is easy to verify directly that it really is a solution, and 
that it involves no complications either at infinity or at the centre of the solid. 
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Introduce a new function W, defined by 

W rr: V — V ( l) 

at all points of space, then we must have 

V 2 W = -Airp (2) 






at all points of space, and, at the boundary, W = and -— — - = 0. 



These last two conditions are equivalent to 

dW = dW = dW 
dx dy dz 



o (^ 



at the boundary, together with one further condition. Equations (3) require that W 
shall have a constant value all over the boundary; the further condition is that this 
constant value shall be zero. 

4. Let F (x 9 y 9 z, X) = be the equation of a family of surfaces obtained by varying 
the parameter X, and such that the boundary of the solid is the surface X = 0. The 
surfaces of this family will divide up the solid into a series of thin shells. There will 
be a contribution from each shell to V t - and also to V . Thus W may be regarded as 
the sum of a number of contributions, one from each shell. 

Let the thicknesses of the separate thin shells be determined by small increments 
in X, say d\ l9 d\ 29 ... , starting from the boundary X = 0. Then we may write 

V,- = Y i (d\ 1 ) + Y i (d\ 2 ) + ..., ........ (4) 

where Y^dX^ represents the contribution to V t - from the shell d\ u and so on. 

Similarly 

V = V (dX 1 ) + V (dX 2 )+ .'... ........ (5> 

Suppose that Y i9 V , and W are being evaluated at a point x\ y\ d on the shell d\ s 
at which the value of X is X'. Then if d\ t is any shell inside the shell d\ S9 the contribu- 
tions to V t - and V from the shell d\ t will be the same; we have Yi(d\ t ) — Y (d\ t ). 
Hence from equations (l), (4), and (5), 

W=V !; -V = {V^xJ-Vo^)} + {YtidxJ-YoidxJ} + ... + {V,.(dx,)-V (dx,)}, (6) 



or expressed as an integral, 

W = I $(V, y\ z\ \)d\. . .. . . . . . . . (7) 



-A' 





5. This form for W satisfies automatically the last of the conditions of § 3, namely, 
that W shall vanish at the boundary. We proceed to determine 3> so as to satisfy the 
remaining conditions which are expressed by equations (2) and (3). 
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Let -r— be used to denote geometrical differentiation with respect too; — i.e., differen- 
ce 

tiation in space keeping y and z constant — and let tt- be used to denote algebraic 

differentiation, i.e., partial differentiation with respect to x keeping y, z, and X constant. 
Then 

dx dx dx 3x' 
and we have 

dw ^ aw ; ^/ aw. . 

dx' dx' dx' lax'j 

- w + £*<**'• *> (9) 

If the point x f , y\ z' is on the boundary we have X' = 0; so that from equation (7), 
-— — = o, and it appears that equations (3) will all be satisfied if 

$(x',y\z',\') = Q, ... ■ (10) 

at all points on the boundary X ; = 0. 

We notice from a comparison of equations (6) and (7) that 

* (x' 9 y\ z\ X') d\ s = V, (d\ s )- V (dX s ), 

and the right hand of course vanishes when x\ y\ z' is on the shell d\. Thus we 
must have, at all points, 

$ (^ */,*', x') = o • • • • • ... ... (11) 

identically, provided that \' has the value appropriate to the point x\ y\ z f . This 
condition of course imposes more restriction on the value of $ than does equation (10). 
Equation (10) was adequate to ensure that the boundary condition (3) should be 
satisfied, but the remark just made shows that for (3) and (2) both to be satisfied — 
i.e., for W to give the true value of V— V , equation (ll) must necessarily be 
satisfied. We shall now assume that equation (ll) is satisfied, and proceed to satisfy 
the remaining condition expressed by equation (2), namely, 

V 2 W= -4tt P . (12) 

In virtue of equation (ll), equation (9) reduces to 

dw = aw 

dx ! dx' ? 
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whence on further differentiation (cf. equation (8)), 

d 2 W 8 2 W d\' 8 2 W 3 2 W dx' 8 , , , , « 

d«' 2 a*' 2 <fa/ ax'Bx' a» /3 ^ a*' lx ' ?J ' z ' h 

so that equation (12) becomes 

_ 4?r = ^ W | ^ W | d " W 
"* p dx' 2 dy n dz' 2 

From equations (8) and (ll), 

3 7 <P (a/, ?/, 2', A') + ^ ^ $ (a/, y', z>, \>) = 0, 



so that equation (13) may be written in either of the equivalent forms : 

1*'™*^ f(d\'\ a , {d\'\ a , (d\'Y\d& /l/A 

/ a$'\ 2 /a$'\ 2 /a$'» 2 

' + 1 ;r-r + 



-4^ = rV 2 «D^-^i ^ \^/ (15) 

Jo 3J^ 

av 

in which V 3 stands for -^— ^ + ~-^ + ^y 2? $ for $ (a/, y\ z f , X) and $' for $ (a?', 2/, 2', A 7 ). 

Thus if $ satisfies either equation (14) or (15), and also equation (ll), then W, as 
given by equation (7), satisfies all the conditions which have been seen to be necessary, 
and will therefore give the true value of V t — - V . # 

* Suppose there are, if possible, two solutions to the same problem, say <£ = $1 and <i> = <£ 2 . Since 
W is determined when the problem is fixed, we must have 

r\> rx' 

Jo • Jo 

so that 

X (^ _ <£ 2 ) dA. = 0, or *i = $ 2 +|f 

, OA 

where 

x (a', y, /, a') - x (&', ?/, %\ 0) = (i) 

for all values of &', y', 5'. Thus, if x is such as to satisfy (i) we may add a term -^ to <£ and still 
obtain a solution of the same problem. A special case in which (i) is satisfied is when 

x (x\ y\ z\ a') = /(a/, y', *', a') {^ (a') - ^ (0)}, 

where / is any function which vanishes identically for the value of A/ appropriate to the point a;', ?/', /, 
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6. As a matter of convenience, involving neither loss nor gain of generality, we 
shall write 

$ Or, y, s, \) = V M/0»> V> *> x )> - . (16) 

in which \Jr(x) * s any function of A. and /(a;, y, 3, X) a quite general function of x 9 y 9 z 9 
and X. Then, from equation (ll), we must have 

f(x 9 y 9 z 9 \) = (17) 

identically at all points, when X has the value appropriate to the point x 9 y 9 z. We 
accordingly have 

so that equations (14) and (15) reduce to 



— 4:7rp = 



/•X' 



, W v^_ <M v){(-:j + (|; + (^}k; • • • <»> 

/a/'Y f^ry fi-H 2 

-4 TP = I ^(x)vy^-f (xO ^- +v vV + ^ (19) 

3X 7 " 

Moreover the family of surfaces (X = cons.) may now be supposed to be determined 
by equation (17), and the boundary will be given by 

f(x 9 y,z 9 0) = (20) 

Thus, to sum up, iff and \j, are such that either equation (18) or (19) is satisfied, 
then the potential of the homogeneous solid of density p whose boundary is 
determined by equation (20), will be given by 



V,-Vn = W= ^ 











xls(\)f(x' 9 y' 9 z' 9 \)d\ 9 (21) 



the value of W being evaluated at the point x\ y\ z\ and X 7 being determined from 
the equation f(x\ y\ V, X 7 ) = 0. 

7. The boundary X = is of course fixed by the solid whose potential is required, 
but we are left with a certain amount of choice as to the disposition of the surfaces 

' Ck' 

nd ^ is any function of X whatever. Replacing ^ (X) - $ (0) by u (A) dX, we find that if # is a 

i j.- I., Jo 



solution then 

'A' 



•♦aWL^M 



will also be a solution of the same problem, 

VOL. CCXV. — A, F 
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X = cons. We shall now limit this amount of freedom by assuming that the region 
at infinity is made to coincide with the surface X = + oo. 

Consider a new function V^ defined, at any point of space x\ y\ z\ by 



.00 



V',= | ^(x)f(x\ y\ z\ X)dX . (22) 





then 



,00 



v 2 v, = ^ oo v y <& 







= -^[Mm)%\ 



A.' =oo 



..... y/iQ J 



by equation (18). Hence, if (as will be the case in all our applications) the value 
of the limit when X f = oo of the term in square brackets is zero, we shall have 



V 2 V y 



^x. 7T « « • . • ■ « . » . • y-Zj^fcy 



At infinity V must vanish, so that at infinity, by equation (21), 



1% 00 



V t . = W = + (X)f (x>, y>, z>, \)d\ = V, 



^0 



Thus Y'i — V t - vanishes at infinity, and satisfies V 2 (Y\ — V t ) = at all points of 
space ; whence (except for a possible singularity at the origin, which will be found 
not to cause trouble) we must have V^ = V t -, so that V' t - is the internal potential. 
Knowing V; and W we find V immediately by equation (l), and have 



V, 



V, 



^{x)f(x\y\z\x)dX (25) 

ir(\)f(a/ 9 y\ /, X)dX (26) 



To recapitulate, the condition that these equations shall give the true values of the 
potential are 

(i) that V; shall be finite at the origin, 

(ii) that \fs (x) ( — ) ~- shall vanish at infinity. 

If these conditions are satisfied, as they will be without trouble in all our 
applications, then equations (25) and (26) will give the potentials. 

8. If f (x, y, z, 0) is the equation of the boundary, the potential at the boundary 
will be 

,00 

is{x)f(x,y,z,X)dX, 



and therefore will contain exactly the same terms in x, y, z as the equation of the 



THE FIGUBES OF EQUILIBRIUM OF ROTATING LIQUID MASSES. 35 

boundary, but with different coefficients. The method is therefore exactly suited for 
the determination of figures of equilibrium of rotating fluid (cf. § 2). 

As a method of determining potentials, the procedure is indirect in the sense that 
we cannot pass by any direct series of processes from the equation of the boundary of 
the solid, as expressed by equation (20), to the general function f (x, y, z y X). We 
must first search for solutions of equations (18) or (19), and then examine what 
problem is solved. 

An obvious case to examine first is that in which f is an algebraic function of the 
second degree. In this case V 2 f is a function of X only, so that the equation for / can 
be satisfied if the last term in equation (18) or (19) is a function of X only. 

Examples of Genekal Theory. 

I. A Sphere. 

9, A quite trivial example may perhaps be taken first, namely that of the sphere 
x 2 + y 2 + z 2 = a 2 . It is seen without trouble that any way of forming the function 
/ (x, y, z, X) will lead to a solution, provided that this function involves x, y, z only 
through x 2 + y 2 -\-z 2 > — i.e., provided the surfaces are taken to be concentric spheres. 
For instance, we may take 

f(x, y, 0, X) = x 2 + y 2 + z 2 -(X-a) 2 , 
then equation (18) reduces to 



A' 





4^ = 6\[r (x) d\ + 2 (X ; -a) yfr (X ; ), 



,1)1' l',/~J* J 

of which the solution is found to be yfr (x) = , H u . 
The potentials are now given by 

1 (X-a) 4 l v / J a 2a dX 3 



V - 

v — 





1 *I£rf\ y +f + z 2 -(\-aY}d\ = M 
(X-a)* L u \ > i T 



? 



where M is written for fTrpa 3 , and r 2 = x 2 + y 2 + z 2 = (X-a) 2 . 



. An Mlipsoid. 



10. It is readily seen that a solution of equation (19) can be obtained by taking 

f( x ,y, t ,X)=-^. + ^L. + ^.-l (27) 

F 2 
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The equation reduces to 



•A' 



i, P = ir (X) 



2 + r^-r + -i \d\ + 4f (X') , . . . . (28) 



a J r\ b + X c +X. 



of which the solution is readily found to he 



Vr(X) 



7T/0 



a&c 



{(a 2 -fX)(6 2 + X)(c 2 + X)} 



i/ * * * ♦ • • • \£J) 



These values are found to satisfy the conditions of § 7, so that the potentials are 
given by 



/»Q© 



, = —Trpabc 






o {(a 2 + X)(6 2 + X)(c 2 H-X) } l/ *V + A ^ 2 +A c 2 + X 



= —nr pabc 



(*CO 



9 « 



^ 



a {(a 2 + X)(fe 2 + X)(c 2 + X)} 1/2 V + X ' 6 2 + X l c 2 ~+A 



1 )d\. 



. ^1 Distorted Mkpsoia, 



11. For an ellipsoid distorted in any way ? and without any limitation (at present) 
as to the distortion being small, assume 

<*>(x, y, z 9 \) = yf, (\) (/+</>) 

where -ft (x) and /have the same meanings as before, being given by equations (27) and 
(29), and $ is any general function of x, y, z, and X. The boundary of the distorted 
ellipsoid is of course 



rv* /J / /y 



-« + -fTJ + ~7> -""" 1 + ( pk ~ — - 0* 

a 2 b 2 e 2 



Equation (19) becomes 



Z.t\2 



— 4 



TTp 



r\ 



o 



O Q O 

-^(x) 'I x + T? + n + ^ 2 ^ (* ^ x + V /l W 



'2x , 30\ 2 , (2y , b>V /2s , dd 
A ray \J3. dy/ \(J d^ 



^2 
4f 



A i -t- B2 C2 3x 



3^> 



in which x\ y\ z\ X ; have been replaced by x, y 9 z 9 X now that there is no danger of 
confusion; and A,"B, C are written for a? + X, 6 2 -f X, c 2 + X. If we further put 



9 9 9 

/y»« 'J/ <^" 

. l_ «/,„. _j_ , 

A 3 B 8 C a 



3x 



ID" 
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the equation becomes 



iirp TtT 



3^ 
3X 



w 



aX/J 



^W{| + | + | + vv}dx 



yfr(\) J4^ 2 + 42~^+ 2 



/a^ 2 



A 3^ 



>#/ 



But from equation (28) 



47Tyo( OT' 



3# 
3X 



or 



DMM+ 



2_ 

c. 



dX + 4 j tu 2 



&)*<*>■ 



so that on subtraction we find as the equation to be satisfied by (p. 



«■ - d * 



-A 



ex, 



*(x)vVA + *(xH«!£ + 4£ + S (gn-a. . (so) 



The equation is too complicated to be attacked directly, but can be effectively 
broken up by assuming a solution 

</> = u+ fv, 

in which u and v are general functions of x 9 y, z, and A, while f is given by 
equation (27). On substituting this value for <f>, equation (30) reduces, after 
considerable simplification, to 



-a' - 



3\/J 



*\{\)V 2 (u+fv)d\ 



+ \{S(\)\4:(1+V) 



X du . du 
Adx 



+ £M S 



X dv dv 
A ox ' d\ 



dcf 



+ ^--i;- s (s + /Sl 



= o. 



and this will be satisfied if we satisfy separately the equations 

^ (x) V 2 (u + fv) dX + 4^ (x) v = 0, 



4(l+t>) 



'y x du du\ Jy, x dy dy 

Adx 3X/ J \ A3x 3x 






0. 



«. • 



(81) 

(32) 



On substituting for /and \js, and writing ABC = A 2 , equation (31) becomes 



4tV 



in which 



= _ f A jv 2 ^ + /V 2 i;+2 

Jo I 



2 ' y( A + B + cJA = 



4x dv 
A dx 



+ 2'W(-r + ^ + 



1 , l\]dx 
U^B + C/Ja 



>IO 



4f 
A 



(— \ - r 

\ A /A «= Jo 



4_ 3i? 
A 3a 



<ZX. 
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Thus equation (31) is equivalent to 



•a r 



1 J \ Adx 3x/J A 



47? 



A=0 



» » * » » [oo I 



This must be satisfied for all values of X, so that we must have (as is clear on 
putting X = in the equation) v — when X = 0. 

Tt will be remembered that the boundary of the distorted figure is given by 



ry& qp> ^2 

a 2 b 2 c 2 9k 



0, 



and it is now clear that A==O reduces to Wx = . Thus the generality of the boundary 
must be involved in the generality of u, and provided u is kept general, we shall 
obtain a general solution of the problem, even if we take the simplest possible value 
for v. The most general way of satisfying equation (33) is to take 



V% + /V^ + 4(2~^ + ^ 

1 A dx 3X/ 



°x 



A^ (34) 



where x ma y be any function of x 9 y, z, and X which vanishes for X or 0, but the 
simplest way of satisfying the equation is to take 



V% + /V=, + 4( 2 || + |) = 0. 



. . • e • t • KdOJ 



In each of these equations the sign of identity ( = ) is used in place of the sign 
of equality, because in the integrand of equation (33) the value of X is not the same 
as the value of X in the upper limit of the integral, which is determined by the values 

of x, y y z. 

12. To shorten the algebra we may change to a new set of variables X, £ ?/, f 
connected with the old variables X, x, y, z by the relations 



X =X, 



x_ 









Differentiation with respect to the new variable \ is given by 



_3_ 



new 



*"" — ~" ~\ jZ*J *" — ™ — * * — — - 

3\ old ' 3x Sec 



3 



X 



X 



where, since x = (a 2 + X) £ we have ^- = f = T , and so 



dx 



A 5 



3 ___ 3 y x 3 

3X new 3X old A OX 



» » # 



. . (36) 
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/ g2 g2 02 \ 

The old symbol V 2 has been taken to mean ( ~~ 4- ^-^ + ^ ) , which in the 



new variables becomes 



1 a 2 1 ^1 1 a2 



12^1^2 > 



A 2 3f B 2 3^ 2 C 2 3f 



and this will be denoted by V 2 ^. 

Equations (32) and (35) in the new co-ordinates, reduce to 

4Mt ,+ 4l) +: 4(t + 4J= (87) 

v>+/vv+4l E0 (38) 

We are assuming that /-I- cf> =■ when X has the value appropriate to the values of 
£ ;?, £, so that in the first equation ^/may be replaced by — <f>> but in the second 
equation this may not be done. 

13. It is convenient, for the purposes of the present paper, to suppose the distortion 
to start from the undistorted ellipsoid, and to proceed in powers of a parameter e. 
Thus we assume 

it = eu 1 + e 2 u 2 + e s u s + . .. , 

v = ev 1 + e 2 v 2 + e*v 8 + ... , 

and in the equation (37) since /+ {u+fv) = 0, it is clear that, when e is small, /will be 



en. 



a small quantity of the smallness of e. As far as <?, equation (37) reduces to -- 1 — 0, 

giving u lm Equation (38) then gives v l ) equation (37) taken as far as e? will then 
give u 2y and (38) will give v 2 ; (37) taken as far as e 3 will give u 3 , (38) will give v z 
and so on. 

As far as e only, equation (37) reduces to -~-~ 0, of which the solution is 
u x = x (i, n> where x is the most general function of £ /;, and £. At the boundary 
^reduces to (u^ = or tox (-1, n, ^) ? so that the generality of the function x enables 

us to deal with the most general small displacement possible. 

At present we shall consider only solutions for which x is algebraic and of degree 
not greater than 3 in g, rj y £. For these solutions equation (38) shows that v x will 
be of degree not greater than 1 in £ r\, £ sa that V\\ = 0, and the equation reduces 
to. 

4 ax ~ ^ " ~~ v a 2 3f B 2 a, a + c 2 ar J Ul 
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or 



&%\ = 






a 2 



a c 



c 



1 
B 



1\ 3f 

,2/32 



1 1\3 S V 



2 / P *'* + I "R &»j 92 + l xC ^ g^^i. 



(39) 



since t; must vanish when X = 0. 



0^1 



Remembering that -^- J = 0, and that f 



eu 



1 + ..., the terms in e 2 in (37) now 



give 



A ou* . di\ 



3X 



3X 



2 



1 /3t^ 



A 2 V3f 



. 1 






-h 



'^ \ 2 



3£/J 



_ lv 1 ° 



o -fcL< 



2 



A 8 3/? 2 






giving on integration 



&Uo = — 



,v 2 



JL 

2 



•A / 



\ 



1 3 2 1 3 2 1 3 2 \ 

Ji 3TL, + g,^ + ^ j ^'dX + u (& * ft 



(40) 



in which « is again the most general function possible of £ /?, f (enabling us to carry on 
the distortion to the second order in any way we please), and the lower limit of the 
integral is taken to be zero simply as a matter of convenience. 

The addition of a perfectly general function « would be equivalent to the superposition 
of a perfectly general distortion (proportional to e 2 ) on to the distortion already under 
consideration. The real object of the present analysis is to be found in its ultimate 
application to the problem of the rotating fluid, and to solve this problem, it will be 
found that w need contain no terms of degree higher than 4 in £ tj, f, this being also 
the degree of the other terms in u 2 . Hence in what follows it will be supposed that 
u 2 contains no terms of degree higher than 4 in £ rj, f. 

A value of v 2 is obtainable from equation (38), but there are, as has been seen, many 
possible forms for v 2i and the most convenient is, in point of fact, obtained by going 
back to equation (34), which in x 9 y, z co-ordinates is 



V\ h+ /V 2 v 2 + A ( 2 f^- 2 + dv * 



A ex 3x 



^\ — 



5x 

3X 



(41) 



where x may be any function of x, y y z, and X which vanishes (to the power of e we are 
now concerned with) both for X and 0. 

Let two new functions w and iv f be introduced, defined by 



t , dw ^ x dw \ 



dx 



V 2 u 



2? 



* a « » 



(42) 



A ,'dw' ^, x div' 
4 (.7: — + 2j 



Jd\ 



A dx 



= -4V 2 w. (43) 
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Clearly since u 2 is of degree 4 in x 9 y, z, w will be of degree 2, so that w f will be 
a function of X only. The term 2 -r- -» — in equation (43) is therefore zero in the 

present instance, but is inserted to maintain symmetry. We now have 



^(l+ s l£)Mt +2 !s)" +4/ (J; +s 



x 3 \ / 



A 3# 
Aw'w 2 — V% 2 — i/V 2 w ? 



so that after simplification, 



V 2 ^ 2 +/ V 2 ( w +/«/) + 4 (^ + 2 £ — ) (« + >'), 



/~~2 



Since / vanishes for the appropriate value of X, ^— ■ will vanish for both X and 

MM 

provided w r is made to vanish when X = 0. Thus ^— will satisfy the condition to be 
satisfied by x in equation (41), and a solution of this equation will be 

v 2 = w-t-fw'. . ■ • • . (44) 

Since v 2 must vanish when X = (§ 11) it appears that both w and v/ must vanish 
separately when X = 0. On transforming (41) and (42) to £ tj, f co-ordinates 
(q/1 equation (36)) and integrating, we obtain as the values of w and w f which vanish 
when X = 0, 



r\ 



w — — i\ V 2 hi;U 2 d\; w r = — -J- V 2 ^ivd\ (45) 

Jo * Jo 

14. Let us introduce a differential operator D, defined by 

d = a i\z + a r)Z + /i_i\z (46 ) 

W AJdf \b 2 B/3^ 2 [c 2 C)d£" K ' 

noticing that, as a function of A, D is purely a multiplier. We have 



aD ^ i a 2 _. 

ax a. 2 af ^ 



and when A = 0, D = 0. 
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The value of v 1 already obtained (equation (39)) is 

and the value of u 2 (equation (40 )) is 

u 2 = -J™ (£ i;, f)-ij g- D^i 2 d\ 

~ i w {£•> *!> 5/"~~"s"^^i • ••••••••••• (48) 



Hence from equations (45) 



•A 
1 I 

w = 



4 V 2 ^ t ^ 2 dX, 
Jo 






'■* ! o dA 



■^(iw-iDw! 2 )^ 



~-^~g JLJCO T""^^. JD Uj_ . . » . . • * . ♦ » . • K~k\j I 



tv f 



= -if V W W ^ 

Jo 

= — i ^r- (— tV D«+^r DX 2 ) d\ 

Jo cX 



- i - TV 3 ,. - l_ "n- 3 ^ 2 



15. This completes the solution as far as the second order of small quantities^ We 
shall not attempt to evaluate u z and v s , as the problems discussed in the present paper 
require a solution as far as e 2 only. 

As far as e 2 , the value of <£ has been seen to be given by 

= u+fv = e(u l +fv 1 ) + e 2 (u 2 +fiv+f 2 w') (50) 

and the potentials can now be found directly from the formula (§11) 

W = f ir(\)(f+<p)d\. 

Jo 

As in § 7, examine a function V' t - defined by 

V',= *(x)(/+#)e*X, 

Jo 

then 

V 2 Y\ = \ ^(x){V 2 f+V 2 (u+fv)}d\ (51) 

Jo 



.00 



Now the value of i^{\) V 2 fd\ is by § 10, equal to — 4tjt>, while from equation (31) 

Jo 
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we have 

*>oo 

yjr{\)V 2 (u+fv)d\= -4:(yfr(\)v) } 

Jo 



00 



Inspection of the values obtained shows that the limit of \fr(\)v when X = co is zero, 
so that equation (51 ) reduces to V 2 V',- = — 4tjt>, and since Y\ is equal to the true 
value of Yi at infinity, and is finite at the origin, Y\ must be the true value of the 
internal potential. Thus the potentials are given by, 



.oo 



V<= ^(A)[/+e(w 1 +A) + e 2 (w 2 +/w+/V)]d\ .... (52) 



V n = 



«oo 

i'(*-)[f+e(ui+f» 1 ) + e 2 (u 2 +fw+f 2 w , )']d\ .... (53) 

*) A. 



in which all the quantities must be transformed into x 9 y> z co-ordinates before 
integration. 

When X = 0, u 2 reduces to Jw (£ j/, f) by equation (40), 



or 



t0 iw &' I' ?) while Wi = x {i > * = * ( J' w ?} 



also v l9 w and w f all vanish when X = 0, so that (cf. equation (50)) 

_ ( x y z \ i 2 ( x y z 

and the boundary of the distorted ellipsoid is 

16. Before proceeding further it will be convenient to examine in detail the 
first order solutions which can be obtained from the foregoing analysis, classifying 
them according to the degree n of the algebraic function u u and, for brevity, omitting 
the continual multiplier e. 

n = 0. Solution is u x = k, v x = 0, <f> = k. 

n = 1. Solution is u x = pg+qr) + r£, v x = 0, <f> =^jr 4- ~ + ^. 

A Jd Kj 

n = 2. u x = af +/3tj 2 +yf + 2frf+2g$g+2hgi, 



v 1 =^i\y% c u 1 dX=^X^^ + 



y 



2n r 



c 2 G 



x% l ^ f\ , q(V 2 l * -f\ , z^ 2 l ^ 



G 2 
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A physical interpretation of these first few solutions can readily be found. For 
the undisturbed ellipsoid of axes ha, kb, he, and origin at x Q , y 0> z , 



3> = 



irpabc 



A 1/2 B 1/2 C 



1/2 



f (x-x ) 2 (y-y a Y (z-z () ) 2 p \ ,, ,. 



and the special ellipsoid which has been under consideration has been that for 
which x = y = z = 0, h = 1. We can change the centre and axes of the ellipsoid 
contemplated in equation (55) by varying x 09 y 09 z Q9 a 2 9 b 2 9 c 2 , and &. If we change 

¥ by an amount §F in equation (55), the change in $ is given by §& — - ^Jrnfrnifr ^ 2 > 

so that y may be regarded as replaced by /-{•</> where $ = — <5P. Thus the solution 
^ = represents a change from F to F—/c; physically it represents a change in 
the scale of the ellipsoid. 

Similarly, if in equation (55) x is changed by Sx 09 y by Sy 09 and z by §z 09 we find 

that $& == -~2^ (X) (^ + ^ + ^\ so that 

</>- z^~ A -t- B -r c 

and the solution n = 1 is seen to represent a motion of the centre of the ellipsoid. 

Similarly, if we put x = y = z = and h = 1 in equation (55) so that $ = ■££ 
and differentiate logarithmically with respect to a 2 , we obtain 

V/a^ 2 "~ 2 a 2 2 A /(a 2 +x) 3 ' 

whence 

Clearly, then, the first three terms in the solution n = 2 represent a distortion of 
the original ellipsoid produced by a change in the lengths of the axes, and it is 
easily seen that the complete solution represents a change of this kind combined 
with a small rotation of the axes. 

n = 3. There are ten terms in the general cubic function of £, r\ 9 £ For the 
present purpose it is convenient to regard this general cubic function as made up 
of a term e^£ and the sum of three expressions such as 

For the solution given by u x = egy£ 9 we have V 2 ^a 1 = 0, so that v x = and 



ABC* 
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It will be shown later that an ellipsoid distorted in this way cannot possibly be a 
figure of equilibrium for a rotating fluid (§ 20). 
For the solution 

^ = £(«f+/V+yf) (56) 

we have 



Vi= -I 



o 



V« Wl ^A - -^g-^ + - 2 ^ B + g-^ 



and the solution is 



as / x* Q y 2 z 2 \ \[x 2 ^y 2 z 2 i\/3aX , /3X , yX\ /K(7 x 

It will be shown later that this distortion leads to the Darwin-Poincare series 
of pear-shaped figures of equilibrium of a rotating fluid. 

n = 4. The analysis of §§13, 14 was confined to the case in which u x was 
supposed of degree not greater than 3 in £ ??, £ But if, in the solution finally 
obtained in § 15, we take u x = 0, which involves taking also v l = 0, we are left with 
a solution (cf. equation (50)) 

<j>^e 2 (u 2 +fw+f 2 iv'), 

in which u 2> w, and w f do not vanish on account of the occurrence of the arbitrary 
function w. And since w has been supposed of the fourth degree, this solution 
gives us the solution of degree n = 4 to the first order, the parameter e 2 replacing 
the former e. Thus the solution of degree 4 is 

u x = a general function of degree 4 in £ ^, f (the old e 2 ^ 2 ) 

This solution is not discussed in detail in the present paper, but is classified here 
with the other solutions for the sake of completeness. An ellipsoid distorted in 
accordance with this solution would give rise to a series of dumb-bell shaped figures, 
which would be figures of equilibrium for a rotating liquid. They would be unstable 
for a homogeneous mass, but the corresponding figures might conceivably be stable 
for a heterogeneous mass (cf. Poincare's remark quoted in § 1 of the present paper). 

17. One point of interest must be mentioned here in connection with the potentials 
derived from these solutions. 

In the potentials arising from the solution of degree n = 2, u± = af 2 , the internal 
or boundary potential will be of the form lx 2 + my 2 + nz 2 , where Z, m, n do not 
involve x, y, z, or X. Since this must be a solution of Laplace's equation, l + m + n 
must vanish, and the potential must be expressible in the form m (y 2 —x 2 ) + n (z 2 —x 2 ). 
All the other potentials may be similarly treated. 
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Making use of this simplification, we arrive at the following scheme for the 
contributions to the internal or boundary potentials of the various solutions up to 
n = 3. Only typical terms are taken ; <p b represents the value of <j> at the boundary, 
Y b represents the contribution from the typical term to the boundary or internal 
potential. 



n = 0. 



fc 



K. 



,oo 



V^ = K \ \jsd\. 



n = 1. 



fc 



tAJ 



a 



2> 



n = 2. (i) <p b = -gp, 



(") <h 



n = 3. (i) fa 



(ii) <t>i 



m) <f> b = 



tAJ 



a 



4> 



xyz 
a 2 b 2 c 2 ' 



xy* 
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2 h V 



tAJ 
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V, 



V re 



tAJ 



o A. 



JUtJ 
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V — — 



^oo 



yj/k \ y 2 —x 2 z 
o 2o£ w A. I Jd 






^ _ /y«^ 



c 



•00 



1 f c&A. 



V, 



V, , = - 



}fr\x_ \x 2 —3y 2 z 2 —y' 



j — 



2b 2 AB\ A ' C 

v 

' \ fr\x \3 y 2 —x 2 3z_ 

2a a A s V B ~C 



2 __ fyfa 

tAJ 



1 r d\. 



-3 kdX. 



18. In any physical application of this method, and in particular in its appli- 
cation to the discussion of rotating masses of liquid, it will be important to know 
what changes are produced by the distortion upon the mass (or density), the 
position of the centre of gravity, and the moments of inertia of the body. These 
changes are given at once by a study of the limiting form of the external potential 
at infinity. 

The potential at infinity of any mass whatever, taken as far as terms of order — 3 , 
has the limiting form 

m + m{x^+y,y + z a z) + 3 (i^ + M^+N« a +2Py« + 2Qza;+2lM 

- — r(L + M + N) +... , 

2r 3 
where 7n is the mass of the whole body, x Q9 y 09 z 0> are the co-ordinates of the centre of 

gravity, and, L, M, N, P, Q, R, are products of inertia defined by px 2 dxdydz = L, 

j j j 

pyz dx dy dz = P, &c. The moment of inertia about the axis of z is \\\ p {% 2 +y 2 ) 
dxdy dz = L + M, and so on. 
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For the solution <p b = k of degree n = 0, the limit at infinity of the contribution 
to the potential is 



§Y=K 



CO 



irpabc -, 2irpabcK 

— (j / \ ___ — . r • • • 

A r 



showing that the distortion involves a change of mass SM. = — 27rpahcK, accompanied 
of course with a change in the inertia terms. 



For the solution <p b = ~ of order n = 1, 

CO 



so that this distortion represents a motion of the centre of gravity by an amount 

tAJQ ~ = ~ "o* CO • 

Solutions of degree n = 2 will clearly involve changes in the moments and products 
of inertia. The limiting potentials are found to be as follows : 

/v.2 2 -j 

(ii) ^ 4 = — , SY X = -%-Trpabc-: -iirpabc 



The first solution does not involve a change in mass, whilst the second does ; both 
distortions affect the inertia. 

For the solutions of degree n = 3, the limiting values are as follows : 

(i) & = J|^> <^ = -f Trpabc^- 

This distortion changes neither mass, centre of gravity, nor inertia. 

.2 „„..2 _ 

05 



(hi) <£ 6 = ~ 5 iV tt = —^Trpabc— —J irpabc 



a *' — • v-r^ >7 «><■#— Sr3 



These two latter distortions move the centre of gravity, but do not affect the 
mass or inertia. 

It is clear, without detailed examination, that the distortions represented by 
solutions of degree n — 4 cannot affect the centre of gravity, but may affect the 
mass and inertia. 
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Figure of Equilibrium of Eotating Masses. 

The Jacohian Ellipsoids. 

19. The condition that a single figure shall be a figure of equilibrium for a 
rotation to about the axis of z is that the centre of gravity shall lie on the axis 
of z 9 and that V & + |-&> 2 (x 2 ~{ y 2 ) shall have a constant value over the boundary X = 0. 
In searching for a series of figures of equilibrium we must add a further condition 
of constancy of mass. 

For the undisturbed ellipsoid, with the notation of § 10, V 6 + i|V (x 2 ±y 2 ) 

^ w} + , {-[ *g * +h ,} + j{- [^ A y (M ) 



x 2 J - ' 7rp 







For this to be constant over the surface, it must be identical with 






where x is a constant. If we put 



2 



2irpahc irpabc 

the equations obtained by comparing coefficients are 



00 d\ 6 

-— : — n = -«> 

AA a 2 



,fl0 dx_ ___ 6_ 

! o AB U b 2 

d\ e 



• <3 O 9 • ft 



. . (59) 



$69800 



'CO 



Jo 



AC c 2 



OOO* • * » 9 



. . . (60) 



• • (61) 



Ellipsoids with Distortions of the First Order. 

20. We proceed to consider which of the distorted ellipsoids can give rise to 
possible figures of equilibrium. 

The solutions of degrees 0, 1, 2 lead to nothing except new ellipsoids, so that the 
inclusion of these distortions could only represent a step along the already known 
series of Jacobian ellipsoids or Maclaurin spheroids. 

Consider next an ellipsoid distorted by the addition of a solution of the type (i) 

of degree 3 (§ 18), say <f> h = e^t^ This distortion, as we have seen (§18), does 

Oj o e 

not affect the total mass or the position of the centre of gravity. The boundary 
of the distorted ellipsoid is 

a 2 b* c 2 V&V 
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while the additional term which has to be inserted on the right of equation (58) is 

f°° irpabc -, 

Hence the additional equation which has to be satisfied, in addition to equations 
(59) to (61), is 

dX e ......... (62) 



Jo 



AABC aW 



Eliminating 6 from this equation and (61) we obtain, as an equation which must 
be satisfied if the distorted ellipsoid is to be a figure of equilibrium, 

AG\ 1 ~(^ + x)(V + X)) = °- ••••••• W 

This obviously cannot be satisfied, for the integrand is positive for all values 
of A. We conclude that the distortion now under consideration cannot possibly give 
rise to a figure of equilibrium. 

21. There remain nine terms for consideration in the general cubic function. 
Inspection will show, or it will soon become apparent as we proceed with the analysis, 
that these fall into three groups, as in § 16, and that the three terms of any one 
group just suffice to give a possible figure of equilibrium when combined with a 
term to restore the centre of gravity to its position on the axis of rotation. We 
shall accordingly consider a distortion in. which the cubic terms are those already 
written down in equation (56). These terms are seen (§18) to move the centre 
of gravity parallel to the axis of x 9 and to correct this we shall add a term (cf. § 16), 

kX , 

~-r- tO U x . 

A 

Thus, for the distorted ellipsoid now under consideration, the boundary will be 

- 2 + fi+- 2 -i+^+^.+r^+VJ- • ■•■••• (64) 

As far as terms in — , the value of the potential at infinity is (cf. § 18) 



r 



-Trpabce- 1| ~ + A |. + A- J + | * j 



> 



so that for the centre of gravity to remain at the origin we must have 



or b 2 & 



K — - -g-l — + 7 2 + 2 ) • . • , W^/ 
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Collecting the terms in Y b as calculated in § 18, we find 
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• (66) 



For brevity in printing, introduce the following notation. Let 
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dX 



A ABC ... 
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•oo 
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A ABC ... 



-ABC * • * 5 



• • # * • * • 



(67) 



so that, for instance, 



J A ■ 
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dX 

AA ? 



I- A A 



-AA 



*>0 



dX 

AA 



And, for the problem immediately in hand, write further 



c i — -Labc — * 



oo 



dX 



A A ~RP ? ° 2 ~~ AAC? ° 3 ~~~ AAB? 



and as before put 



CO 
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9 then equation (66) becomes 



[v.+ v («*+»•)] 



irpabc 
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{J A -n) + y 2 (J B -n)+z 2 J G 



OLX r 2 
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The distorted ellipsoid will be a possible figure of equilibrium, as far as the 
first order of small quantities, provided the right-hand member of equation (68) is 
identical with 



irp 



ahce i + ¥ + $- 1+ex ( a $ + v^ +r ^ + ^\- • • • (69) 



Equating coefficients, we obtain 



J A -» = -» J B -w = p Jc = - 2 ' (70) 

^(c 2 + c 3 )-2pC 3 -^- 2C2 = 0^, (71) 



Iaa + tIsIab+^Iao +^a = 0~, (74) 



3ot 



and from these equations, together with equation (65), we must eliminate or obtain 
the coefficients. 
If we put 

?=«' h* W (75) 

then equations (71) to (73) reduce to 

*'(cj + c 8 )-0'c8—/c 8 = — a! t (76) 

CO 

90 

a'(-Sc s )+/3'(3c 3 + c i )-y , c 1 = ^- 2 /3', (77) 

CO o 

90 

a'(-3c 2 )-p'c 1 + y'(3c 2 +c l )= 4^ 2 y', (78) 

CO 

whilst on substitution for k from equation (65) and J A from (70), equation (74) 

T*Pnl]OP^ to 

K(3a , I A A+iS'lAB + yiAc) = *wa s (3a , + /3 f + y) (79) 

We have now to deal with four equations (76), (77), (78), and (79), and have 
to examine whether, and how, these can be satisfied by values of a\ (3\ and y. 
Since there must, from general physical principles, be an equation of some sort for 
points of bifurcation (whether capable of being satisfied by real values or not), we 
are led to suspect that these four equations are not really independent. 

H 2 



D Zi 
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Equation (79) was obtained by the elimination of k from two equations (65) and 
(74), each of which expressed in effect the condition that the centre of gravity 
of the mass should be at the origin; in fact, equation (65) was only a short way 
of arriving at the value of a-, which would in any case have been given by equa- 
tion (74). We therefore expect that equation (79), derived from (65) and (74), 
will prove only to be an identity of which the truth is involved in the three other 
equations (76) to (78). And, as a matter of procedure which is entirely at our 
choice, we shall elect first to solve equations (76) to (78), and then to verify the 
truth of (79). 

The elimination of a\ /3', y from equations (76), (77), and (78) gives a determinant 
which on expansion reduces to 



a 



l h V L-% + p + -^ ) (ci c 2 + C! c 8 + 3c 2 c 3 ) 



20 

a 2 



{^(tf + c 2 ) + c 2 (3a 2 + c 2 ) +c 3 (3a 2 +& 2 )} + (~^f = 0, (80) 



and this is accordingly the equation giving points of bifurcation on the Jacobian 
series of ellipsoids. 

22. Two identities of importance are the following : 

~~ 3 /I 



J A + Jt5 + J 



BT^C 



From equation (70) 



»00 



\A B C A 



+ ^ + ^ ) d X 



o 3x \A, 



.00 



AA \A b c 



2 X 



dx 



3 / 1 



abc 



> • • (81) 



3X \AA 



dX=2J A . . (82) 



J 



A 



n + 



e 

a' 



2? 



(83) 



T , 6 

It ^^K^^mM f if I I „.. I, ,n 

eJ B — ■ n -+- 2 , . 



(84) 



J 



8 



c 



c 



ii 



• » « » 



so that on addition, by the use of (81), 



, , . . . • KoOj 



2 



abc 



Xfi 



1 1 _, 1 
2 + 6 s 



2 n + i — -r Ti i" « ; > 



(86) 



giving 6, and hence J A , J B , and J c , in terms of a, b, c, and n. We further have 



,00 



Xati — 



I 



AC 



\d\ 
AAB 



a Ja~—£ Jc 

2 2 ? 

a —c 



,oo 



a 



o l(a 2 -6 2 )AA (a 2 -6 2 )AB. 



■c?X 



a 2 J A ~& 2 J 



B 



a 



3 -6 2 



d a 



• • W 



. . (88) 
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and I AA is given by equation (80) as soon as the values of I A b and I AC are known. 
Substituting for J A , J B , Jc ? from equations (83) to (85), we obtain 

T _ T - a2 T _ 1 I ° 2 2e \ /oM 

1 A B ~ U 9 1 A C — " 2 2 n > J-AA — "~ 3" i — 2 1 71 21* ' * V° y / 

a — c [ a —& a \ 



Finally we have 



t'i - — 



A ABO c 2 -W 



1 1 \ A d\ 1 /j j \ 



C 2 



— ~2 2(IaA~Iac) ? C 3 — 2 7 2 (IaB — Iaa). 

c —a a —b 



With this material it ought to be possible to find the points of bifurcation from 
equation (80). As, however, Darwin's results are available, it will be sufficient to 
make use of his results and merely verify that his quantities satisfy equation (80), as 
they are in point of fact found to do. 

23. Darwin's values, calculated for an ellipsoid such that abc = 1, are 

a = 1*885827, b = 0*814975, c = 0*650659, 

n = -^~ = 01419990, 
whence, by equations (86) and (89), 

%e = 1 ~ U 1 = 0'2068037 
a 2 b 2 c 2 

I AB = 01419990, I AC = 01611871, 1 AA = 0-0711382, 

c x = 0-07967602, c 2 = 0-02874219, c 3 = 0-2450100. 

With the use of these values, equations (76), (77) and (78) become 

, 01216184a , +0-02450100/3 / + , 02874219y / = (90) 

0-07350300a' + 0-1970290 / 8' +0-07967602y' = 0. .... (91) 
0-0862266a' +0-07967602 ( 8' + 0-3835217y / =0 (92) 

The values of a', ft, y, are, of course, indeterminate to within a common multiplier. 
The simplest set of values, obtained by cross multiplication of the coefficients of 
equations (90) and (91), is 

a' = - 0'003710945, & = O^OOH 43630 y' = Q-000595338. 
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If we substitute these values in equation (92), we find 

0'0862266a / + 0'07967602/3 / + , 3835217y / 

= -0 , 00031998 + 0"00009111 + , 00022829 = -0*0000005. 

The fact that the error occurs only in the seventh place of decimals adequately 
verifies Darwin's calculations, but the tendency of small errors to accumulate in 
computation is forcibly illustrated by the circumstance that in the above equation 
the final error is as much as one-six-hundredth part of the whole value of the middle 
term. 

With the values just obtained for a', /3', y\ I find 

la 2 . (3a / + j8 / + y / ) n = -0-00094878, 
\a 2 (3a 7 Iaa + ^Iab + tIac) = -0'00094885, 

verifying that equation (79) is satisfied, again as far as the sixth place of decimals. 

24. With a view to subsequent computations, it is convenient to take a standard 
set of values such that a! = —1. These values are found to be 

J a! = -1, $ = 0*3081810, y' = 0*1604294, 

and with these we have, by equation (65), 

K = - \ : (Sa! + p + y') = 0*506278. 

These numerical values substituted in equation (64) will give the equation of 
PoiN cake's pear-shaped figure as far as small terms of the first order. 

The Pear-shaped Figure Calculated to the Second Order. 

25. The question as to whether the pear-shaped figure is stable depends upon the 
change effected by the distortion upon the angular momentum of the ellipsoid. But 
(c/*. § 18) the first-order distortion so far considered can be easily seen to produce no 
effect at all upon the angular momentum of the figure. It is therefore necessary to 
proceed to terms of a higher order, and we now consider terms of the second order. 

The first-order terms have been found to be given by 

u x = £ (af + fy 2 + yt? + ic), ........ (93) 

with (cf. equation (47)) v x = —^Du^ The value of u 2 will be given by equation (48), 
in which u x is to be assigned the value (93), and w will be taken to be given by 

It has to be shown that this value for w makes it possible for the figure distorted 
to the second order in this way to be a figure of equilibrium. 
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It will be noticed that with the value we have assumed for w, the value of u 2 
becomes a function of £ */, £ of degree 4 involving only even powers of £, */, £ ; its 
value is 



Au 2 = 2?(a£ a + fy*l yt 2 + k)[ 3a( 1 - ^) + #( jjj- - l) + y(l - i 



2 



+ L£ 4 + M* 4 + N? 4 + 2Zi; 3 f 3 + 2m . . (95) 

The values of j% and of jf V are easily seen to be similar in form, so that all the 
second-order terms in <p are of this form (of. equation (50)). Multiplying by \fs(\)d\ 
and integrating from to oo, we obtain as the terms of the second order in the 
potential an expression of the form 

— wpabce 2 (c n # 4 + c 22 y 4 + c 33 z 4 + c 12 x 2 y 2 + c 2 $fz 2 + c 31 z 2 x 2 + d x x 2 + d 2 y 2 + d z z 2 + d 4 ). 

If this figure can be a figure of equilibrium at all, it will be for a rotation differing 
only by a second-order quantity from that of the original ellipsoid. Let us suppose 



2 



that for it 7- = n + e 2 n" ; then at the boundary, as far as e\ 

-[T ft +K(^+ifl] 

= Trpa&c { x 2 ( J A — n — eV) + y 2 {J B —n— e 2 ?i") + z 2 J c • 

+ e 2 (c n # 4 4- c 22 ?/ 4 + c 33 2; 4 + c 12 ^V 2 + c 23 g/ V + c zl z 2 x 2 + c£ 3 «x 2 -f d 2 y 2 + c£ 3 2 2 + d 4 

+ third-degree terms in e, the same as before} (96) 

At the boundary, 

so that for the figure to be one of equilibrium, the right-hand member of equation (96) 
must be identical with 

! /vj2 j .2 ~2 / sy,2 2 2 

a 8 & 8 c 8 6 4 c 4 cW a 4 6 4 



» 2 qp ^ 

-.+2qf-+2r- t 



+ 2p^ i +2q y u +2r Z - i +s)\. . (98) 
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Upon equating coefficients of the terms of degree 1 and 3 in x 9 y, z 9 we obtain 
exactly the same equations as were obtained before when the terms in e 2 were 
omitted. Thus no alteration is produced in the cubic terms until the solution is 
carried as far as e 3 . 

The equations obtained by comparison of the coefficients of x 2 , y 2 , and z 2 are : — ■ 

J A - n - eV + e 2 d x = 0[- a + ip-X 

and two similar equations. In these the terms in e 2 and the terms independent of e 
must be equal separately. The latter terms again give equations (70), while the 
terms in e 2 give 

7 // 1 OP / Q Q \ 

\A)-\ '"" IV — "o T • * * • • • • • • • \i/w) 

2 a 4 v } 

(1/2 "**** IV —~ 2~ ■j-Z . • • . • • • « • . (-LUUJ 

3 ' — ' "2" 4 * " * * ' ' * ' * * \ / 

Finally, upon equating terms of the fourth degree, we obtain the six equations : 



c n 



= i-Tg-; C a9 = i-Tr; c 33^i--r. • . (102) to (104) 



a c 



c 2S—2^^'y C 31~2^4' c i2~2^ 4 - • . (105) to (107) 

We have seen that equations (70) to (73) must still be true, so that a, b, c, a, /3, y, k 
will be the same as before. 

26. At present there are eleven quantities to be determined or eliminated, namely, 
L, M, N, Z, m, n 9 p, g, r; s> and n". The equations giving these quantities can be 
simplified in the following way, 

By an argument already used in § 17, it appears that the terms in e 2 on the right 
hand of equation (96) must be harmonic. Thus we must have 

6c n 4- c 12 + c 13 = 0. (108) 

6c 22 4- c 23 4- e 2l = 0. .......... (109) 

^^33 ' ^31 ' ^32 ^^ V ^IIOj 

d x 4- d 2 + d 3 — (llL) 

These are, of course, not new equations to be satisfied ; they reduce to identities 
when the calculated values are inserted for o u , c 12 , &c. 
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With the help of (ill), equations (99) to (101) may be replaced by 

w " = -^(ij+$ + £)' •••••••• (112) 

d 1 -dE, = 46(^-i), (113) 

U/g - ==: 2^(7 — ^ , . . . . . . . . . . . .(114) 



while from equations (102) to (110) we obtain 

3L m n , v 

3M + I + ^.-n hl6^ 

* c 4 a 4 

3N + I+^ =0 (117) 

C 4 & 4 <^ 4 \- LJ -'/ 

* 

These three equations with equations (102) to (104), namely, 

On = ±^> ( 118 ) 

°22 ~ 4 7 8 J • * \ ±J -^/ 

c 33 = i -?, (120) 

may be used to replace the group (102) to (107). 

27. We proceed to evaluate these quantities in detail. The values of %i x and v x from 
equation (47) are 

U X = a£* + ptf + yg* + K& (121) 

Vl= ~* iX \^ + m + ^k]> ••••••• ( 122 ) 

and the value of u 2 is already given by equation (95). 

For convenience in computation we shall combine all the terms in u 2 which are 
independent of X in the first two lines, with the similar terms in the last line ; we 
accordingly write 

- 1 (3af + M +y ?+ K y -pr» 2 - ^ set 

+ L£ 4 + M^ 4 + N£ 4 + 2 ^f + 2mf f + 2nfn + 2 (pf + qn + H 2 ) + s 

= L'£ 4 + MY + N V + 21' J? + 2m' ?e + 2n'? n * + 2 (p'f + q'f + r'f ) + s' 

~~~ CO \ £*, W, >/• • » • • • • • » • * » • • • « • • * • \ 1 juO / 

VOL. CCXV. A. I 
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The value of u 2 is now given by 

+ l(Saf + ^ + y^ + K y + ^ + ^£l + o,'(i, n , Q ■ • (124) 
whence, by equation (45), 



Aw = 



2 



^ LvA u 2 dX 



o 

2^2 / c.2 2 \ / ji2 o£2 

+ P \ '2 -n2 +0 A TD MY I 2 ,-12 ^ 



a A , .^ ^ B3 . - AB y , ^ C2 . AG 

+ ^(|| + ^) + «r( A § + |>/3r(|3 +A ^ + i B X 

+ K 1 | 2 + K^ 2 +K 3 <; 2 +K 4 , (125) 

where the K's are constants, to be chosen so as to make w vanish when X = 0. 

o o o 

The value of V\ , n f (4w) is derived from that of Aw by replacing £ 2 , j? 2 , f 2 by -r%, ■=-., ™ 

1 " " A B O 

and omitting all other terms. Tims we have, again from equation (45), 



W = - i V\ „ f (4w) c£x 



_ jL J j r, a 3/3 3y _^ 3ot/3 3 ay /3y 

- 4 | "2 A8 + 2 AB 2 + g AC s "» A2 g + A 2 C + ABC 

+ 2 f — - 2 + 2-^, + - A + T^ + "p 2 " 1 " 5 |' • • • • (126) 

in which K 5 is a new constant, chosen so as to make v/ vanish when X = 0. 
On collecting terms, we obtain 



where 



4(w+/w / ) = P 1 £ 8 + P a i, 3 + Psf + P 4 (127) 

8 A 2+ 8 B 2+ 8 C 2+ *AB + 4 AC + 4 BC 

■ oil' ,3 M'A 8 N'A ! Z'A x m' lJL w' 
+ t3s A + «"T?~ + " 8 "~C 5 " +4 BC ,l¥ "c + l4 B 

+ liK 1+ i^+i^ + |K 6 A (128) 
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p _ i7 a2 B , o3 j3L , 3 r!B „ a a/3 3 oyB lJL j3y_ 

^ 2 ~ 8 A 1 + 8 AB + 8 AC 2 + 4 A 2 + 4 A 2 C + 4 AC 

' L'B ,, S M' 8 N'B +ll Z' im'B , llW ' 

+¥ A F + 3 * B +¥ "C^ +l¥ C +¥ AC~ +l¥ A 

+i^ + liK 2 + i^+iK 5 B (129) 

p _ l7 a 2 C g /3 2 C „g y 2 3 a/3C q3 ay , -■ i /3y 
■fs - l-« A3 + 8 ^ +^8 ^ +4 ^ig +^4 JT+U £g 

+i^+i^ s + iiK,+iK,c (130) 



r * "U 2 AB + AC/ + A + B + C + 4 



1^-5^1 ,3 _^!_ ,3 y 2 | 3a/3 3ay_ _/3y 

2 A 3 2 AB 2 + 2 AC 2 A 2 B + A 2 C ABC, 



« +2 ^ + | + | + | + K 6 ) (181) 



2> 



If the value of <f> is taken to be efa + e 2 <f> 2> we have as the value of <p 
4<p 2 = 4w 2 + 4/ (w +/W) 

+ fe U B + C/ + *V A + A + A/ + A 

+ L'f + MY + N'£ 4 + 2 W + 2mTi a + 2» W + 2 {p'g 2 + q'n 2 + r'£ 2 ) + s' 

+(Af+B, 2 +cr-i)(p 1 r+B !! , 2 +p 3 r+p 4 ), (132) 

and we have already supposed (§25) that 

2^-— = c n x* + c 22 y* 4- c 3S z* + c 12 x 2 y 2 + c iz x 2 z 2 + c 23 ?/ V + d } x 2 + c^?/ 2 4 d d z 2 -t d 4 . (133) 
A 

I 2 
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Hence, by comparison, we obtain as the values of the coefficients 



c 
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+ 
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3 






P 

JL 9 
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^ J A 2 « ■ -M. j- 2y \ 2^ P4 _ Pj_ 
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»ao 



A \AB S ~W B ~ B 5 



^ /2/cy 2r' P 4 _ P/ 
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00 



A \A 
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140 



141 
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143 



Evaluation of Certain Integrals. 

28. It is clear that before these coefficients can be evaluated, certain integrals 
must be calculated of the types J B ca..., Iabc.» where the notation is that of §21 
(equation (67)). 

The values of J A , Jb, Jc have already been evaluated in § 23, as also of I A a, Iab, 
I AC and I AAB , I AAC , Iabo- We also have 



2 



T _ If j c 2 T _ nb 

1bc - y^? B J/-c 2 ° ~ V-~<? 



s = 0'3916228. 



These 10 integrals will form an adequate basis from which to calculate all the 
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others. The required integrals can all be obtained by successive applications of 
formula? of the following types, which can be verified without trouble : 

(a 2 - b 2 ) J AB = J B - J A ; (a 2 - b 2 ) I AB = I B - I A , &c. 



31 



BB 



Zo n — J-nn — J- 



B 



BC X AB > 



a 2 J 



AA 



2 



(2n+ 1) J A n + l ™ <% J )C ^ a 2 n Ja»B 



J a — Iaa, OCC. 



AA? 



J A n c ' 



(2U+ l) I A n + l = 2 J A n — I A n B — I 



A n O 



A great number of the integrals can be calculated in two or more ways, and owing 
to this circumstance it is possible to provide very complete checks on the accuracy of 
calculation. Two complications are worthy of mention. 

In the first place one must consider the ordinary cumulative error of all prolonged 

computation. Since b 2 is nearly equal to c 2 , any error present will be increased when 

f(l> 2 ) — f( c 2 ) 
we evaluate a new quantity of the type J2 2 — I consequently where a quantity 

can be evaluated in two ways, the one in which division by b 2 —c 2 is not involved has 
been taken to be the true value, and the one derived by division by b 2 —c 2 , has been used 
merely as a check, and has generally been found to differ in the sixth or seventh place 
(or near the end of the computations even in the fifth place) from the other values. 

Secondly, if the 10 integrals used as base were known with perfect accuracy, the 
checks ought to be satisfied fully except for the error in the last one or two figures. 
But, as has been indicated in § 23, the 10 integrals are not themselves perfectly self- 
consistent, so that different methods of computation will lead to a difference of the 
final results comparable with the errors in the basic integrals. 

The following table gives the values I have selected as the best for the various 
integrals required. I have not thought it necessary to record the checks or estimate 
the probable errors here, as a much more searching test of the accuracy of the whole 
computation can be provided at a later stage. 



J A = 0*2583003, 
J AA = 0'05262769 ; 
Jbb = 0'6516017, 



J = 1-8401326. 

J B = 07647290, 
J AB = 0-1751040, 
J B0 = 0-8813026, 



Jo = 0-9769708, 
J AC = 0-2293883, 
Joe = 1-2044842, 



Jaaa = 0-011873224, J AAB = 0^04234772, J AAC = 0-05641920, 
Jabb = 0*1647550, J ABC = 0-2254075, J ACC = 0'3112352, 



Jbbb = 0-6830283, J B bc = 0*9537991, J BCC = 1'9011148, J C oc = I '9011148. 
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J 
J 
J 
J 
J 



Jaaaa = 0-00281568, 


Jaaab : 


= 0-01053694, 


Jaaac — 


0-01421838, 


Jaabb = 0-04232384, 


Jaabo : 


= 0-05842981, 


T _ 

"AACC — 


0-08133328, 


Jabbb = 0-1791995, 


Jabbc : 


= 0-2518508, 


Jabco = 


0-3563870, 


Jaccc = 0*5074646, 


Jbbbb 


= 0-788915, 


"BBBC = 


1-124336, 


Jbboc = 1-611794, 


<J BCCC : 


—~ Zi ozi/ yOj 


"CCCC = 


3-361221, 


aaaaa = 0-0Q06881964, 


» AAAAB 


= 0-002669728, 


J AAAAC 


= 0-003639559, 


aaabb = 0-01099071, 


J AAABC 


= 0-01528665, 


« AAAOO 


= 0-02142202. 


aabbb = 0-04732648, 


•J AABBC 


= 0-06687764, 


« AABCO 


= 0-09532252. 


AAcoc = 0-1360144, 


J ABBBB 


= 0-2108169, 


« ABBBC 


= 0-3016727, 


abbcc = 0-4340725, 


« ABCCC 


= 0-6273273, 


"ACCCC 


= 0-910874, 


I A = 0-9215282, 


Ib = 


= 1-3322118, 


Io= 1'4265252, 


I AA = 0-0711382, 


J-ab - 


= 0-1419990, 


Iac = 0-1611871, 


Ibb = 0-3319454, 


Ibc = 


= 0-3916228, 


Ice = 0-4670439, 



Iaaa = 0-01040244, I AAB = '0"02450100 : 

Iabb = 0-06567633, I ABC = 0-Q7967602 : 



I AAC = 0-02874219, 
Iacc = 0-09762467, 



Ibbb = 0-19794510, I BBC = 0*2478016, I BC c = Q-3131750, I CC c = 0*3996337, 



-l-AAAA = 

J-aabb = 
Iabbb = 
Iaccc = 

J-BBCC = 



0-001859707, 

0-01423688, 

0-04573358, 

0-0963965, 

0-2714534, 



-AAAB — 



-AABC — 



-ABBC — 



-BBBB — 



-BCCC — 



0-004874751, 


Iaaac 


= 0-005853759, 


0-01761100, 


Iaacc 


= 0*02198620, 


0-05813153, 


J-ABCC 


= 0-07452920, 


0-1590430, 


Ibbbc 


= 0-2070218, 


0-3590070, 


1-cccc 


= 0-4811801, 



Evaluation of the Coefficients c ll9 c m .... 

29. It will be noticed that the coefficients c n , c 12? ... are linear in a, a/3, ... , 
I/, M 7 , N', ... , p\ q\ r\ s' so that the various contributions may be calculated 
separately and independently. 



Contributions from Terms in p' 9 q', r\ s\ 

As regards the contributions from these coefficients, we may take (cf equations 
(125) and (126)) 



K x = K 2 = K 3 - ; K 4 - - 2^ K 5 - 

a • 
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whence we obtain as the contributions to P 1; P 2 , P 3 , P 4 {cf. equations (128) to (131), 

4 A + B C 4 \ali 6 a B VC/' 

It now appears that p\ q\ r\ s ! contribute nothing to the values of c n , c 12 , . . . , c 23 
{cf equations (134) to (139), p. 60.) Their contributions to 4c? l9 4d! 2 , and 4d! 3 are as 
follows : 

4^ = 



o A \A 2 A 



'°° c?X /2p ; p r X g ; X r' X \ 
.o A \A 2 a 2 A 2 6 2 AB c 2 AC/ 



so that the contributions an 



&d 1 — 2p J AA niAA - 77-1-AB ^J-ACj 

or b 2 c 2 



t „r „j 



Ad - 9n f T - ^ T - # T r T 

4a 2 — zg Obb 5 -Lab tI-Lbb £-Lbc 



a" o 3 



c 

4ta s — Z/7 JcC 5 -LAC*"" T^J-BC 2 -LcC« 



a 2 ^ b 2 nKJ & 

Since this part of the potential should be harmonic, we ought to have d 1 + d 2 + d s = 
(d/l equation (ill)). This is clearly the case, in virtue of the identity 

2a J A a — -Laa + J-ab + -Lac- 

I have verified that these identities are satisfied by the values in the table opposite, 
and the contributions are found to be 



j ~t ~j 



4:d x = 0'303186^--0*141999^--0161187 ~ 

a b c 

U 2 = -0'141999-^ + 0'533622^ 2 -0'391623x 

ar b 2 c 2 ' 

4d 3 = -0-161187 £ 9 -0 , 391623? 5 +0 , 552810- 9 

a 2 b 2 c 2 

Contributions from Terms in U , M', Nj, V, m', n'. 
30. As regards these terms, we may take {cf. equations (125) and (126)), 

1_ \a 2 b 2 c 2 

J ~ W V + c 2 

K 4 = 0, 

T / 7/ TZ T / 7/ 

a be a ar b c 
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and the contributions to P 1? P 2 , P 3 , P 4> are found to be (cf. equations (128) to (131)). 

/ T ' \ 2 V \ 2 

p __ _ 1 [ H y ±± A__ 1 y ° A 

4 ~ 4 V a a* A a &V BC. 



3IA m'X n'X 

\ a 2 A c 2 G b 2 B, 

I'X 3M'\ n'X N 

,c 2 G b 2 B a 2 A, 

' I'X , m'X 3N'\N 

b 2 B VA c 2 G 



p _ _ ZdJjA mA 71A\ _,p 
X 1 — 1 ..2» T J2C\- ' 7.2TJ ] - tLL 1 



P - _ LL + ™L* + !L± _BP 



P 3= _(±^ + ^ + ^>)_CP 4 . 



We find as the contribution to d5 I? (equation (140), p. 60) 

1 ~J AU AV~Jo A I A + A 2 W 2 A + c 2 C + & 2 Byj' 

from which it is easily verified that d l + d 2 + d i = 0, as it ought to be. 

In virtue of this relation it is only necessary to evaluate two of the contributions 
4d v id 2 , 4d 3 , but I have calculated all three directly from the table on p. 61, so as to 
obtain a check on the amount of error involved from the causes mentioned in § 28, as 
well as a check on the accuracy of my own computations. The values I find are 



L' _„„*„„ M' n . nnnnnn N' 



Ad, = 0-85375 ~ - 0-073783 ~ - 0-089893 , 

■ 4 (r c 4 



0-054134 J-, + 0-072732 4^ + 0-Q59701 n 



be c a a o 



T f M 7 NT 7 

id 2 = -0-041149— + 0-244072 4=- -0194266—, 

a 4 6 4 c 4 



7/ / / 

4- 0'051969 T |~- 2 - 0*054134 ^ - 0*005568 -^ , 

6V cV a 2 6 2 



T 7 M 7 ~K f 

4c? 3 - -0-044227 ±L - 0*170277 ~ + 0*284157 —, 

a o c 

a f f 

+ 0*003091 ~ - 0*018600 ^ - 0'054134 — 

b 2 c 2 cV a 2 b 2 

From these figures the value of d r fd 2 + cZ 3 would be given by 

4(^ + ^2 + ^)= -0'000001±L + 0*000012 ^_ - 0*000002 ^- t 

a b c l 

+ 0*000026 ~ - 0*000002 ^ - 0*000001 n 



b 2 c 2 a 2 c 2 a 2 b 2 



This check gives an idea of the amount of error involved. It illustrates the tendency 
of the errors to accumulate in terms where b's and c s are plentiful, and to be absent 
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where as are present. This is a consequence of the method of procedure explained 
in §28. 

In calculating the contributions to c u , c 12 , . .. , I have (unwisely) taken the variables 

1/ M' 

a" b 4 

contributions need have been calculated, but I have calculated five out of the six quite 
independently, so as to have two independent checks on the computations. 
For the contribution to c n , I find 



to be I/, M 7 , ... , instead of — , -r- , ... . In virtue of relations (108) to (110), only three 



.00 



— I ^ / _L_ 1 

11 ~ Jo "Ml 1 A 8 

= h (J AAAA "~ 77~2 IaAAA + 7T~Z I AAA ) + M ( ~y^ IaAB "~" ^T2 -^AABB ) 

+ W \8c* Iaac ~ 8? Iaaoc J + V \4&V Iaao ~ 4? Iaabc ) 

+ m ' (icV Iaaa ~ 4a 3 1aaac ~ 7 Iaaac ) + n> ww Iaab_ i& Iaaab 

The other five contributions can be similarly written down. It then can be verified 
algebraically that the three relations (108) to (l 10) are satisfied ; this provides a check 
on the method and the algebra, but it hardly seems worth exhibiting it here, as the 
numerical checks to be given later provide simultaneously checks on the method, the 
algebra, and the computations. By independent computations I find for the contribu- 
tions to c n , c 22 , egg, c 12 , c 13 , the following : 

4c u = 0'001359233L' + - 01278935M / + 0-04066156N / 

+ - 01515472r-0-01251121m / -0 , 006581149n / 
4c 22 = , 000446082L / + , 1490178M / + - 2y ? 8016lN / 

-0 , 390933lZ / + 0-00681885m / -0 , 00912350w / 
4c 33 = , 000576222L' + 01129554M / + 0-385800N / 

-0-3972137/ / -0-0l767165m' + 0-00509336w' 
4c 12 = ~0 , 003687295L / -01465540M / + , 2013742N / 

-0-0642884r-0-03593784m' + 0'06239405» / 

4c 13 = -0'004468110L' + 0-0698178M / -0-4453435N / 

-0-02663960/' + 0-11100470m / -0-02290710» / . 
From these 

4(6c n + c 12 + c K! ) = -0 , 000000007L / -0 - 000000lM / + , 0000000N / 

-0 - 0000003Z / -0"0000004m' + , 00000000w', 

I have not calculated c 23 independently, but a check is provided by the comparison of 
the values of c 23 deduced from the above figures by the use of equations (109) and (110) 
vol. ccxv. — A. K 
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separately. For these values I find the following, the upper value being that derived 
from equation (109) : 

4c 23 = 0*00101080 0747552 1*86947 

I/— M 7 — N' 

0*00101078 0*747550 1*86946 

2*40989 0*0049753 0*0076531 

+ V — m r — n. 

2*40992 0*0049754 0*0076515 



Contribution from terms in a 2 , a/3, . . . , &c, 

31. For the computation of these contributions it is convenient to take the standard 
set of value obtained in § 24, namely 



a' = 4 = -1, P = fi = 0*3081810, y = ^ = 0*1604294, 



a 2 *' " & 2 "> / ^ 



* = -X^a' + ^' + y ) = 0*506278. 
These give the values 

a = -3*556343, /3 = 0*2046890, y = 0*06791892, 



a 



2._io«/2/fi7Kh70 /02 A.AjTon^^ 2 



12-647573, /3 J = 0-0418976, y J = 0*004612981, 



a /3 = -0-7279442, ay = -0-2415430. /3y = 0-01390225. 

Then, as regards the terms in a 2 , a/3, ..., &c., including ;c, we find (cf. equations 
(125) and (126)) 



^! 4. 3 @L , 3 r! , 6o^ 6oy /3y 



*>—<&+¥+&)— rw+v+rt 



2 o /O \ ^2 



\a c ar a b / a 






TT — _ _ Mir _i_ ti i I — 



TT -„|L5ii 3 J^L , 3 j/L l §a£ , §oty /3y K x K 2 K 3 

R ~~ l 2 ~ a 6 2 a 2 6 4 2 aV a 4 6 2 aV aW a 2 6 2 c 2 / 



''TZ* TZ" "17" \ 

a o c 



giving on substitution of numerical values, as regards terms in a 2 , a/3, . . . only 

K x - -19*91885, K 2 = 0*1102214, K 3 = 0*04221664, 

K 4 - 0*3603665, K 5 - 3*556845. 
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The values of P 1; P 2) P 8 , P 4 originating from these terms are 

A 3 - A 3i^ 4 8 Jl +1 8 £2 +J-8 £* • t> 4 ^ +^4 AC +U BC 



A 2 ! 4 A + 4 B + 4 € 



5 — 1 fi7^ |o3 fi 2 , 3 y 2 1 o3 gfi , ,3 <*>y , ii ffyl 

B 3 ~ AW ^8 A 2 +^8 B 2 +8 C2 +^4 AB +4 AC + 14 BC r 

+ B 2 l 4 A + 4 B +4 C 4 5 

P« _ 1 [1.7. °* 2 t 3, /3 2 . o3 y , 3 ^ - 03, ^y ■ a Py\ 
C 3 ~~ AC 2 1 8 A 2 +8 B 2 + 8 C 2 +4 AB ~^ 4 AC + 4 BCJ 

1 J 1 K x ! K 2 1 jl K 3 1 -r^ 



p _ J^a j8 y \ ,jr _l/i5 ^l , 3/3 3y 3a8 3ay __/3y 

4 ~ \A 2 AB AC/ 4 4 \ 2 A 3 "*" 2AB 2 + 2AC 2 A 2 B " h A 2 C + ABC, 

4 U + B + c +K 7 

1 have evaluated independently the three contributions to d u d 2y d% and find for 
them 

4<i x = 15a/cJ AAA + 3 / 8/cJ AAB + 3y/cJ AAC + K 4 J A 

— (26|a 2 J AAAA + 2i/3 2 J AABB + 2iy 2 J AACC + 7^a/3J AAAB + 7^ay J AAAC + lj/3y J AAB c) 

— ( 1 2"Ki J AA + 2^2 J AB + 2" 1^3 J AC + 2^5 ^ A/? 

and there are corresponding values for the contributions to id 2 and 4<i 3 . The 
numerical values found by independent computations are 

4dj = 0*032398, id 2 = -0-011175, 4d 8 - -0*021223. 

As a check on the computations, the sum of these contributions ought to be zero ; 

we find 

4 (di + da + dg) = 0-000000. 

The six contributions to c ll9 c 129 ... have also all been calculated separately. The 
values obtained are 

4c u = 0*0722711, 4c 22 - 0*0271162, 4c 33 = 0*0274693, 
4c 12 = -0*215756, 4c 18 = -0*217869, 4c 23 = 0*053056, 

as regards the checks which ought to be satisfied by these values, I find 

6c u + c 13 + c 13 = 0*000002, 6c 22 + c 12 +c 23 = -0*000003, 6%+^ + % = +0*000003. 

K 2 
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Formation and Solution of the Equation. 

32. We proceed to form and solve the final equations. It is convenient to deal 
with the equations in c ll9 c 12 ... first. 

By comparison of coefficients in equation (123), we obtain 

= L' + 50*01206. ..... (144) 

- M' + 0*01178109 .... (145) 



L 


T / , 9ot 

— JLi f 2 

a 


- 10 KCL 


M 


— J\lL 4" — t> 

a 




N. 


2 

— i_lj -J™ 2 

a 





= N'+ 0-001297113 .... (146) 
l = y+ftt =V + 0-003909145 .... (147) 



2 

a 



m = m'+-^ + % -5/cy = m'- 0'3538936 . ...» (148) 



2 



n = n ' + i^ + M _5^ = n' - r0060520 . .... (149) 
a b 

p^y + l^ c -5/c 2 = #'- 2*800420 ..... (150) 

a 

q = q'+^ = a' + 0-02913935 .... (151) 

-Mm •JL, J~\t "^ 

a 

r= ,/ + r5 = r > + 0*009668880 .... (152) 

a 2 v 



s = *'+- 2 = /+ 0*07207330 .... (153) 

a 2 K ; 



These coefficients can be checked by comparing the sums of the coefficients in the 
two values of «' (equation (123)) ; i.e., taking f 2 = ^ = f 2 = 1. For the difference of 
the sums which ought to vanish, I find 41*86192 — 41*86191 = 0*00001. Using the 
values just obtained, 

~ L = (T002585675I/ + 01293149, 

a? 

r, M = 2125367M 7 + 0*02503914 > 

6 y 

■^ N = 12-87541^ + 0*01670087. 

Collecting the various contributions to c n , c 22 , c 33 , equations (102) to (104) now 
become 

0*0025856751/ + 01293149 

= 0*0013592331/ + 0*01278935M / + 0*04066156N / 

+ 0*01515472^-0*0125112 W~0*006581149?i / + 0*0722711 3 
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2 , I.25367M' + 0'02503914 

= - 000446082L / + - I490I78M / + 0-2780I6lN / 

-0*390933lZ' + 0*00681885m'-0*00912350ft' + 0*0271162, 

12*8754lN' + 0*01670087 

= 0*000576222L' + 0*1129554M' + 0'385800N' 

-0*3972137Z'-0*01767165m' + 0*00509336ra' + 0*0274693, 

while equations (115) to (117) give 

1/ = -23*52190m'-9*55672»'-32*07326, 

M' = - 0"8204313^-0'01162649to / -0-00329142, 

N' = - 0'1354301/ / -0'00472373m / ~0-000154830. 

The solution of this set of six equation is 

I/= -6172711, M'= -0-01761613, N' = -0*002105706, 

V =_ 0-006053622, m' = +0*5865530, to' = +1*659250 . : (154) 

No check is needed on the accuracy of these values beyond the fact that they 
satisfy the equations. On substituting directly into the equations, I .find that they 
are all satisfied accurately to the last place of decimals. 

The corresponding values of L, M, N, I, m, n (see opposite page) are 

L= -1171505. M = -0*00583504, N = -0*000808592, 

Z=_ 0*002144477, m= 0*2326594, n= 0*653198,. . (155) 

and the values of c u , c 12 , ..., calculated directly from equations (102) to (107) are 

4c u = -0'03029132, 4c 22 = -0*01240160, 4c 33 = -0*010410956, 

2c 2S = -0*01121809, 2c ;!1 = 0*04245102, 2c 13 = 0*04842301. (156) 

These ought to satisfy the checks afforded by equations (108) to (110). In point 
of fact, I find 

6c u + c 12 +c 13 = +0*00000007, ' 6c 22 + c 23 +c 12 = +0*00000012, 

6c 33 + c J3 + c 23 = +0*00000006. 

33. The first use which must be made of these numbers is to determine the 
contributions to d 1} d 2 , d s , evaluated on p. 64. We have by separate computation, 

0*085375 ^-0*073783-fr-0*089893--... =-0*3412367, 
a b c 

0*041149^+... = + 0*1672652, 



a 

T / 
Ju 

a 



0-044227^+... = + 0-1739738, 
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Hence collecting all contributions to 4c£ 1? 4d 2 > 4c£ 3? we obtain as their total values 

4d x = 0-303186^0-0141999 2.-0161187 ™0 - 308839, 

a o c 



id 2 ^ -0141999^ + 0*533622 2.-0-391623 ^ + 0*156090, 

a be 



.t 



4c? 3 = -0*161187^-0-391623 2, + o*552810 ^+0*152751 

a 



■ 2 -— J, 2 v — c 2 



From the values already obtained for p, q y 'r, we can transform equations (112) to 
(114) into the following : — 

^ + r. + ^) -01163013^+0*6227300 2- 2 + 0'9769708 --0*0419475, 

a 4 i> 4 cV a 2 h 2 c 2 

dd^d^ = 26^-^) =0'2326027^-l'2454600 |-2-0'2378029 5 

4d 3 = 20™ =1*9539416 ^ + 0*04462528. 

3 c 4 c 2 

On substituting the values of 4c£ l5 4d 2 , and 4c? 3 , these reduce to 

0*212582 £- 2 + 0'569839 f 2 + 0'230436 \ - 0*227126, 
a 2 h 2 c 2 

0*161187 -^,+ 0*391623 ^+1*401132 ~ = 0*108126, 
a b 2 <f 

0*116301 £,+ 0-622730 2-+o*976971 ~ = 0*0419475-4^". 
a 2 b 2 c 2 

The solution of these equations is found to be 



a? 



= l'428257 + 32'04689w // (157) 



<£ = -0 , 111620-11797935tc // (158) 

- 2 = -0' , 0559390-0'3891163n // (159) 

c 

and this solution has been verified by direct insertion into all the equations. 

34. This completes the solution of all the equations, and the determination of all 
the coefficients except s'. 
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On collecting the values of K x , K 2 , K 3 , K 4 , K 6 , from pp. 62, 63, and 66, and inserting 
the numerical values already obtained, we find 



K x = -(^-+a + -i) -19-91885 = 28'26823, 



K 3 =-(^ + ^r+y +0-1102214 = -0-2624723, 



K 3 = - ™+f- 2 +^fR 0-0422 1664 = -Q-0986790, 



31/ n' m /N 
a 2 b 2 c' 

V , 3M' , V 

„a 2 b 2 & 

W + ^ + 3N A 
a 2 b 2 c' 

K 4 = -2—1LJL +0-3603665 = -0-900332- 19-85984»", 
a b c 

K 5 = f 2 I {+S T f i + 3-556845= -2770998. 

To evaluate s f we have to examine the form assumed by the potential at infinity. 
The additional terms in the potential, as far as terms in — , produced by the distortion 

are readily found to be 



—irpabc -(s f — fK 4 + x 3 5K 5 ) 5 
r 

and if Sm is the additional mass produced by the distortion, this must be identical 

(STY) 

with — . Hence, if s r is determined by the condition of constancy of mass, we have 

s f = |K 4 ~ T 3 r) K 5 = -0*230755-13*239893 n" 
giving {of. equation (153)), 

.8= -0*1586814-13"239893^. 

Discussion of the Figure. 

35. The boundary of the pear-shaped figure, as far as the second order of small 
quantities has been found to be 

O O O / f> o o 

a b c \ a b a o ac a 



+ ie 



iJLx* , Mg/ 4 , Nz 4 , 2ly 2 z 2 , 2mz 2 x 2 , 2nx 2 y 2 , 2px 2 , 2qy 2 , 2rz" , \ ' , * 

j-e I — - — V —77, — h— 7-+ 74 , -i rr 1 ttt- + — -7- + 77- H t- + 5 = . (loO) 

\ a 8 6 8 c 8 &V c V a 4 & 4 a 4 6 4 c 4 / ■ ' 



In this equation all the coefficients have been determined ; the coefficients p, q $ r and 

« 
2 

s have been found to involve n", defined in § 25 by the equation -^- = n + e 2 n", while 

27T/0 

the remainder are pure numbers. 
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Let us put 



P ~ Po +Pi n 



tt 



and similarly for q, r, and s, and let us put e 2 n" = f ; then the equation may be put in 
the form 

2 2 ^2 2 2 ^2 ^ 

~2 + f2 + ~l~ 1 -+^ a— 4-/3 -4n + y-!-^ + -o 
a 2 6 2 c 2 \ a b a 2 & 4 ' aV a 2 y 

, ! 9 /La 4 , My" , N* 4 , 2ZyV , 2mzV , 2na?Y , 2p () ,x 2 , 2g ?/ 2 2r 2 2 

' 4 \ ^8 ' 7,8 ~ ^8 ' 7,4^4 ~ „4^4 ' ^,4/vl ' _4 T 7,4 ' Jk ^ *<- 



a 8 ¥ c 8 6V c 4 a 4 a 4 6 4 a 4 6" c 



4 

2 o_.2 o_ ^2 



+ xs" ( — r~ ~! — T^t — I 1 — h §i ) — . . . . . . (161) 

For any values whatever of e and f, provided only that they are sufficiently small, 
equation (161) will give a figure of equilibrium. If we put e = 0, but retain £ the 
equation becomes 

which is an ellipsoid of semi-axes a\ b\ c\ given by 

a ' 2 - 1 4- if f'-Si-L-l-Q ) — - l J-Af/^4- A« \ — - 1 4-+f^- JL +s ^ 

1 " 2S l 9 ' 2*1 /> 72 ~~ l ' 2M72 1 2*1 h 2 ~~ L ' 2M j ' 2*1/? 



a 2 ^ W. A \r V ^ \6 2 z V c 2 ^\c 



or, numerically, 

a n IV 2 r f2 

%;= 1 + 1271347£ t~ 9 - 1~9'20894<;, --- - l-3T)0453f. 

It is at once clear that as f varies this ellipsoid coincides with the various Jacobian 
ellipsoids near to the standard ellipsoid. 

If we put f = but retain e in equation (161) we obtain equation (160) with 
n f/ — ; i.e. we obtain a series of figures of equilibrium all having the same angular 
velocity as the standard Jacobian ellipsoid from which they are derived. 

The two series of figures obtained by putting e = and f = in equation (161) 
may be represented by the two intersecting straight lines POP 7 , QOQ ; in fig. 1, the 
point of bifurcation being of course represented by the point O. The general figure of 
equilibrium represented in equation (161) is, however, arrived at by assigning values 
to both e and £ these values being limited by the condition only that e and f shall be 
so small that e 3 and f s/a shall be negligible. Thus the figures of equilibrium given by 
equation (161) are represented by all the points inside a certain rectangle ABCD 
surrounding the point 0, They do not fall into two linear series, as it seems to be 
tacitly assumed by Darwin and Poincare that they will do. 

36. The circumstance that the two linear series lose their identity and become merged 
indistlnguishably into a general area seems to be predicted as a direct consequence of 
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Poincare's general analysis, coupled with the linearity of the equations (V 2 W = —4^, 
&c.) which lead to figures of equilibrium. For the vanishing of the Hessian (A == in 
the notation of Poincabe # ) which expresses the condition that a point of bifurcation 
should exist, expresses also the condition that the two linear series should be merged 
with an area of linear series as they approach the point of bifurcation, t 

There is, of course, no question that in the neighbourhood of the point O two linear 
series do actually exist, such as may be represented by the lines POP', ROE/ in 
fig. 1 ; this is abundantly proved by Poincare's general argument. % What is now 
maintained is that an expansion as far as e 2 only, does not suffice to reveal the 



Jacobian. Ellipsoids 

(unsbable) 



R 



P 1 




•^ 



Jacobian Ellipsoids 
(stable) 



Fig. 1. 



directions in which these linear series start out from the point of bifurcation. So 
long as pur vision is limited to the interior of the rectangle ABCD in fig. 1, we can 
know nothing of the direction in which the line OR starts out from O. And the 
whole difficulty is merely one introduced by the artificial method of expansion in 
powers of a parameter ; as soon as this artificial method is abandoned the rectangle 
ABCD shrinks to an infinitesimal size, and the curves POP 7 and ROR ; become merely 
two lines intersecting in the point without any complications. An exactly 

* " Sur Pequilibre (Tune masse fluide animee (Tun mouvement de rotation," ' Acta Math.,' VIL, p. 259. 
t Cf. footnote to p. 74. 
I Loc. cit.y § 2. 

VOL. OOXV. — A. L 
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analogous situation arises in considering the directions in which lines of force start out 
from a point of equilibrium in an electrostatic fields 

37. A more interesting illustration of the difficulty will be found in an investigation 
of the figures of equilibrium of rotating liquid cylinders which I published in 1902.t 
In this paper the equation of the cross-section of a figure of equilibrium corresponding 
to a rotation « is supposed to be expressed in the formj 

^-(l~^-)(lfo + ^ + ^ 2 + &0 3 +...) 5 (162) 

where £ 0? £, £ 2 , £ 3 .... are functions of ;? ; £ rj are complex co-ordinates given by 
£ = x + iy, q = x—iy, and x, y are ordinary Cartesian co-ordinates measured from the 
axis of rotation. The quantity 6 is a parameter, analogous to the e of the present 
paper, measuring distance from the point of bifurcation at which the pear-shaped figure 



2 
0) 



and the elliptic cylinder coalesce. At this point of bifurcation, 1 — - — = f , so that 
£ = £ 0J of which the value is shown to be 

£ = h±W(h 2 -~Wa>), 
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If Y is the potential of an electrostatic field, the equation of a line of force will be 
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where I, m, n are direction-cosines. Two lines of force will meet in a point of equilibrium (just as two 
linear series meet in a point of bifurcation), and the condition for this is 

a v __ ev _ ay _ 

da; 03/ c# 

Let £c , ^o> ^0 he a point of equilibrium satisfying (ii.), then, if e is a small quantity of the first order, the 

point 

x + Xe, y + fie, z + ve 

will, as e varies from zero upwards, trace out a line passing through % y , z . The condition that this 
shall be a line of force is, as far as first powers of e, 

dv ' dv dv 

dx dy Sz 

X ~~ fJb ~~ v ' 

and this is satisfied (analytically) because of equations (ii.) for all values of A, /x, v, . Thus, as far as first 
powers of e, there are as many lines of force through the point of equilibrium as there are values of 
X, /x, v ; an infinite number. But on going as far as e 2 , it becomes clear that there are only two true lines 
of force through this point. The condition that a point of equilibrium shall exist is also the condition 
that, if the approximations are not carried far enough, there shall be the confusion of an infinite number 
of lines appearing to satisfy the condition for a line of force, and the analogous statement is true for 
points of bifurcation and linear series of figures of equilibrium. 

t " On the Equilibrium of Eotating Liquid Cylinders," ' Phil Trans.,' A, 200, p. 67. 

I Lot. cit, equation (71), p. 86, 
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The calculation of £ presents no difficulty, and equation (162) as far as £ only is 
shown to be the equation of a pear-shaped figure. On calculating £ 2 its value is found 
to be of the form {cf. § 22 of the " Cylinders " paper), 

180 3 2825 4375 _ 8 , , . /48 , 1984 _ 3 



6= T"^""^'""2r , " 8+ --' + * , fe 9+ 
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where S 2 is analogous to the n!' of the present paper ; to be exact the rotation for any 
value of 6 is supposed given by 



1-:^- = % + $»& + **&'+.... 



Again, then, as far as 6 2 there is a doubly-infinite series of figures of equilibrium, not 
two singly-infinite series. 

In this earlier and simpler investigation, it was an easy matter to carry the 
computations to the third, fourth, and fifth orders of small quantities. It was found 
that the equations giving £ 3 for a figure of equilibrium could not be satisfied so long 
as $ 2 was kept indeterminate, they could only be satisfied for one special value of S 2 , 

namely $ 2 = . After having determined the value of S 2 in this way, it was 

possible to investigate the stability, and the pear-shaped cylinder was found in point 
of fact to be stable. What is important in connection with the present paper is that 
it was not possible to determine the stability of the pear until after its figure had been 
determined to the third order of small quantities. # 

38. The work of Poinoare can hardly be compared in detail with the investiga- 
tion of the present paper, because Poinoare tacitly assumes the whole point at 
issue ; namely, that it is possible to determine the beginning of the pear-shaped series 
from an investigation of figures of equilibrium going only as far as second-order 
terms. The work of Darwin admits of detailed comparison, because Darwin's 
work claims actually to have effected the determination which I am compelled to 
believe, after my investigation, to be impossible. 

It will be clear that any extra condition in addition to the conditions that the 
figure should be one of equilibrium will provide an additional equation which will 
reduce the doubly-infinite series inside the rectangle ABCD down to a singly-infinite 
series represented by a straight line. For instance, the condition that the angular 
velocity shall remain constant would reduce the rectangle to the straight line QOQ' ; 
the condition that the angular momentum should remain constant would reduce it to 

* The previous investigation on cylinders and the present one on three-dimensional bodies follow- 
widely different methods ; the present paper is in no sense a translation of the former from two into three 
dimensions. The two papers were written at an interval of twelve years, and I had hardly referred to the 
former paper in writing the present one until after I had encountered the difficulty of not being able to 
determine the stability from the second-order figure ; I then discovered that precisely the same situation 
had arisen in my former investigation. 

L 2 
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some curve at present undetermined. Hence, if the present paper is sound, we 
should anticipate that Darwin obtained a linear series instead of a rectangle of 
configurations by the introduction of some adventitious condition not necessary to 
equilibrium. 

39. Darwin supposes his pear defined as far as the second order by the equation 

00 

t = c-eS 3 -2// S/, (163) 

1 

where t gives a measure of the surface displacement, e is a parameter analogous to 
our e (although not numerically the same), and the quantities fl are coefficients which 
must vary as we pass along the linear series, but are constants, as is also e, for any 
single figure of equilibrium. The quantities S 3 , S/ are ellipsoidal harmonics. Darwin 
supposes e to be a small quantity of the first order, and the coefficients f£ to be small 
quantities of the second order. The energy E of the distorted ellipsoid will differ 
from that of the original ellipsoid from which the displacement t is measured by a 
small quantity SE which will involve e, f s and Sa> 2 , where §w 2 is the increase in the 
value of co 2 for the distorted figure. 

The first order displacement of the ellipsoid will be represented by the first terms 

T = c~~ eS 3 , (164) 

of equation (163), and this is supposed to be a possible figure of equilibrium with 
Sco 2 = 0. The corresponding value of §E will be of the form 

SE = ae 2 9 (165) 

and the coefficient a must therefore vanish if the original ellipsoid was one 
corresponding to a point of bifurcation. 

If we do not at present assume that a = 0, the value of §E arising from the second - 
order displacement (163) will be of the form 

SE = ae 2 + (3e* + y e 2 f+6f 2 +^ (a + be 2 +fc) . .... (166) 

in which a single term in f has been taken as being typical of all the terms in the 
coefficients//. The condition that the displacement (163) together with an increase 
§w 2 in co 2 shall give rise to a figure of equilibrium is that E shall be stationary for all 
variations of e and f ; it is expressed by the equations 

A (£E) - 0, A (SE) - 0, (/ = //, &c.). .... (167) 

Expression (166) is the same inform as that given by Darwin (' Coll. Works/ 
vol. 4, p. 349), except that he omits the term ae 2 from SE ; and equations (167) are 
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the same in form as those from which Darwin obtains the conditions of equilibrium. 
The equation «- (SE) = 0, written out in full, becomes 

a + 2/3e 2 + yf+b -^- = 0. (168) 

2irp 

If second-order terms are neglected equation (168) reduces to 

and if a. is taken equal to zero, equation (168) reduces to 

2/3e*+ 7 f+b~=0 (170) 

Darwin's method is in effect to replace equation (168) by equations (169) and (170). 
This, in my opinion, introduces one limitation too many on the values of the variables, 
and so reduces the doubly-infinite series of figures of the second order to two singly- 
infinite series. Equation (169) must undoubtedly be true if second-order terms, are 
neglected, but we may take 

a = X, 2/3e 2 + yf+b^=-- - X (171) 
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Trp 



where X is a quantity of the second order entirely at our disposal, and still satisfy all 
the conditions necessary for equilibrium. I think it will be found that Darwin's 
procedure in effect introduces the spurious condition X = 0. Darwin's equations are 
of course sufficient to ensure equilibrium ; what is maintained is that they are not 
necessary and so do not disclose all possible figures of equilibrium. 

In this way I believe it will be found that Darwin has limited himself to one linear 
series of figures (X — 0) instead of the doubly-infinite series represented inside the 
rectangle ABCD in fig. 1. If this series should happen to run on continuously at 
the edge of the rectangle with the true series ROB/ in fig. 1, then of course Darwin's 
investigation would stand. But no reason suggests itself, and certainly Darwin 
(not foreseeing the complication of the doubly-infinite series) gives no reason, why this 
should be the case. For some value of X the two series will run on continuously, but 
there is no assignable reason why this value should be X — 0. 

What, then, will happen if we try to carry Darwin's approximation on to third- 
order terms by his method ? I think it will be found that his series comes to a dead 
end before the third-order terms are arrived at, precisely as if it ran on to the 
boundary of the rectangle ABCD in fig. 1, and could get no further. If the displace- 
ment goes as far as third-order terms, £E* must go as far as sixth-order terms, and 
will contain terms of orders 2, 4, and 6, say 

oE = SE 2 + (5E 4 + aE 6 . 
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The equations necessary for equilibrium are 

^(SEi+SEt+SEi) = 0, |.(£E 2 +ffi 4 +£E 6 ) = 0, (/ = //', &c), (172) 

while the equations provided by Darwin's methods would be 



k-(£E 2 ) = 

ce 


= o, 


g>((iE 2 ) - 


= o, 


|w = 


= o, 


g>(^E 4 ) = 


= 0, 


|« - 


= 0, 


"Z^AoI-JqJ — 


= 0, 



and it will readily be seen that there are more equations than can be satisfied by the 
variables at our disposal. The conclusions of this section are put forward with the 
utmost diffidence, but to the present writer they seem inevitable. 

40. Assuming that no glaring error has been made, the present investigation seems 
to indicate that so far the knowledge we have as to the stability of the pear-shaped 
figure amounts to absolutely nothing. The required knowledge can only be obtained 
by carrying the figure of the pear to a still higher approximation. In the parallel 
investigation on C3dinclers it was found that the stability could be examined as soon 
as the figure was determined as far as third-order terms, and doubtless the same will 
prove to be the case in the present problem. Fortunately the method of the present 
paper is one which lends itself to indefinite extension, limited only by labour of 
computation, so that it ought to be possible to proceed to third-order terms and 
determine the stability of the pear, and if the pear then proves to be stable, to pro- 
ceed to higher orders and so examine the series of pear-shaped forms. In the previous 
investigation on cylindrical figures it was found that an expansion as far as fifth-order 
terms gave a good approximation to the pear-shaped figure up to a stage where it 
was obviously just about to divide into two separate masses. It is only in the hope 
that I shall be able to carry the present investigation further, that I have ventured 
to put forward the present somewhat lengthy piece of work which has so far led 
only to such negative and disappointing conclusions. 
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